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ABSTRACT 


The Physics Department at the Naval Postgraduate School is developing a 
concept to overcome the problems that keep present rail guns from being practical 
weapons. The rails must be replaced often if the rail gun operation is to be continuous. 
Replacing the rails in present rail gun configurations is time consuming. The Physics 
Department’s design concept uses a rectangular barrel as part of the solution to the 
problem of replacing the rails. The projectile will require flat surfaces to maintain 
electrical contact with the flat rails and aerodynamic stabilization because of the lack of 
angular momentum. This thesis develops one possible model of a projectile for a 
rectangular barreled rail gun, which could be used to replace the standard five-inch gun 
found on most warships. The proposed projectile is successfully modeled as a five inch 
projectile with flat areas planed onto opposite sides and long chord, short span fins 
attached in a cruciform configuration. The computer programs used to develop die 
projectile model are included to allow evaluation of alternate configurations. 
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L INTRODUCTION 


This thesis models an aerodynamically stable projectile for a rectangular barreled 
rail gun. The three computer programs used in modeling this projectile are provided so 
that the models can be updated to meet changing requirements. 

The first section of Chapters II and III present the reasoning and assumptions for 
the selections made from the available alternatives regarding that chapter's subject. The 
second section of Chapters II and IE develop the equations used in the computer model 
associated with that chapter. Chapter IV uses the assumptions from the previous two 
chapters and begins with four sections of equation development. The final section of 
each chapter analyzes the output of the computer program developed in the previous 
sections from that chapter. 

T his thesis is multidisciplinary, with advisors in both the Aeronautical 
Engineering Department and in the Physics Department. Readers from different 
disciplines are reminded that they may find more detail relating to their particular subject 
area than they feel necessary, but readers from other disciplines may need that detail to 
understand what is being done. 

A. BACKGROUND 

A rail gun uses the interaction between electrical currents and magnetic fields to 
accelerate a projectile. The rails in a rail gun are electrically charged with the projectile 
completing the electrical circuit between them An electric current flows through the 
r ails and projectile creating a magnetic field. The magnetic field interacts with the 
flowing electrons, which accelerates the projectile out of the barrel. The acceleration is 
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increased as the magnetic field is increased. The magnetic field increase as the electric 
current increases. 

The Physics Department at the Naval Postgraduate School has a rail quick change 
concept for a rail gun. The large electrical current used in rail guns to achieve the large 
accelerations to attain large muzzle velocities in present rail gun designs also vaporizes 
some of the surface of the projectile. The vaporized material settles on the rails, fouling 
the barrel in about five firings or less. The rails must be removed for cleaning before the 
rail gun can be fired again. Present rail gun designs take thirty minutes of longer to 
remove and replace the rails, assuming clean replacement rails are available. 
Conventionally fired guns can maintain a rate of fire of several shots per minute for 
extended periods. A thirty-minute cease-fire between five shot salvoes is not an 
acceptable rate of fire for a rail gun weapon system. 

The Naval Postgraduate School Physics Department’s rail gun design concept, 
among other things, will shorten the time required to replace the rails and thus increase 
the firing rate in comparison with existing rail gun designs. The computer models 
developed in this thesis provide graphs for the analysis of the aerodynamic effects 
imposed upon a projectile by this alternate design as physical dimensions are changed. 
The resultant projectile performance can also be analyzed. Only those details of the rail 
gun design effecting the design of the projectile, such as barrel shape, will be discussed 
here. Rail gun design is such a broad field of study that all other details of the concept, 
such as power requirements and building techniques. Are left to other investigators. 

The most important feature of the new rail gun concept is its rectangular barrel. 
The projectile, in common with all rail gun projectiles, requires an electrical contact 
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surface capable of withstanding the immense electrical currents used in rail guns. The 
magnitude of these electrical currents depends on the weight of the projectile, the desired 
muzzle velocity, and the length of the barrel. The electrical contact surfaces of the new 
design's rectangular barrel will be flat rails on opposite sides of the barrel. Therefore, to 
maintain good electrical contact the electrical contact surface on the projectile will have 
to be flat. The rectangular barrel, having no rifling, will not impart the rotation that 
provides the angular momentum that usually keeps projectiles stable in flight. The lack 
of angular momentum adds the requirement of some form of aerodynamic stabilization 
for the projectile. 

B. BASIC THEORY 

A variety of aerodynamic theories, presented by S. S. Chin [Ref. 1] and by 
William Pitts, Jack Nielsen, and George Kaattari [Ref. 2], were used in development of 
the computer models used herein to compare projectile shapes and related performance. 
Slender body and nose/body theories from Reference (1) were used to determine the skin 
friction drag, form drag, boattail drag, and base drag of the body of the projectile. Airfoil 
theory, also from Reference (1), was used for analyzing small fins added to the projectile. 
Figure (1) on the next page shows these stabilizing fins attached to the projectile body. 
These fins compensate for the lack of spin stabilization normally provided by the angular 
momentum caused by the rotation of the projectile. The fins, however, also add drag. 
Wing/body interference factors from Reference (2) were used in the equations that 
combine the aerodynamic characteristics of the projectile body and fins to provide a total 
drag coefficient for the equations of motion and the aerodynamic coefficients for the 
stability study of the complete projectile. 
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Figure (1) Projectile Schematic 

3.6 2.4 

H-in ->|e--18.1 in-—►Kin-H ^.34 





Numerical integration theory is required for the computer model that calculates 
projectile trajectory and impact energy. As a projectile travels through space, many 
factors interact. At any point in time the velocity is changing due to the effects of drag, 
which in turn changes the distance the projectile will have traveled at the next point in 
time and the velocity it has when it gets there. To further complicate matters, the drag 
itself changes with the velocity and with the air density, which is changing with the 
vertical position of the projectile. With so many interactions between variables a closed 
solution to the integration of the equations of motion is difficult, consequently numerical 
integration schemes are required to solve the equations of motion. The Euler-Cromer 
method was selected for the numerical integration of the equations of motion because of 
its stability in long calculations. The Euler-Cromer method, occasionally referred to as 
the last point method, uses a loop, shown in simplified form in Figure (2) as a flow chart. 
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The program loop takes the final position of the projectile from the previous time step 
and calculates the velocity at that point and then uses that velocity to calculate the next 
point in the projectile's trajectory. The accuracy of the integration for most numerical 
integration methods is increased as the time steps are made smaller, but the round-off 
error often increases with each calculation. The Euler-Cromer method is stable in long 
calculations because the round-off errors tend to cancel each other so the calculation does 
not diverge. 

Figure (2) Euler Cromer Numerical Integration Flow Chart 
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II. SIZING 


There are many different ways of fitting a projectile in a rectangular barrel. At 
the time of this writing the design of the rail gun is still very conceptual, thus there are 
few physical limits placed upon the projectile. 

A MATLAB computer program, provided in Appendix (A), is used to model the 
different possible projectile nose/body configurations as the radius of the projectile is 
changed. This MATLAB computer program will be referred to as the configuration 
comparison program. Since some of the design alternatives chosen may become invalid 
as the rail gun design is formalized, the equations for the configuration comparison 
program have been generalized. Future users of the program will be able to change 
inputs, such as the nose fineness ratio or the explosive payload of the projectile model, to 
suit new requirements and the model will still be correct. 

A. CONFIGURATION ALTERNATIVES 

The configuration choices made from among the different projectile design 
options, e.g., ogive verses spherical nose, body diameter, etc... will determine the 
projectile's physical dimensions. The physical dimensions, such as cross sectional area, 
will have a direct effect on the aerodynamic characteristics of the projectile. 

1. Electrical Contact Surface 

The projectile requires flat electrical contact surfaces while in the mil gun barrel. 
Two methods to do this are immediately evident; place the projectile in a flat-sided sabot 
or put flat sides on the projectile. Flat sides on the projectile for the electrical contact 
area are used in the model because of the two major problems with sabots outlined below. 


7 



A sabot is a spacer placed around a projectile to allow the projectile to be fired 
from a larger or differently shaped barrel. Sabots are usually light and often used in 
conventional guns when using a larger powder charge than would otherwise be possible 
behind a smaller projectile to increase the projectile's muzzle velocity. A sabot seems 
attractive from an aerodynamics consideration because the projectile can remain truly 
cylindrical. A cylindrical projectile should have reduced airflow separation and surface 
flow field turbulence. The sabot in the rail gun, however, has to conduct the electrical 
current required to accelerate the projectile without melting. T his will require a 
substantial piece of metal, which equates to a large mass and thus a large inertia. T his 
large mass will have to be accelerated with the projectile and the large inertia means it 
will take time to stop once it has left the barrel. 

Accelerating the sabot becomes an energy consideration. The energy at the target 
is half the impact mass times the impact velocity squared. Since a sabot separates from 
the projectile immediately after leaving the barrel its mass is not present at projectile 
impact so the energy used to accelerate the sabot is wasted. 

The large inertia associated with the sabot is a problem because the projectile will 
separate from the sabot as they both exit the barrel of the rail gun and the sabot will land 
somewhere down range. The light plastic sabots in the Phalanx gun, for example, have 
comparatively low velocity, high drag, and low inertia so they travel only a short distance 
and fall harmlessly. The sabot from the proposed rail gun's projectile will have hi gh 
drag, but it will have high velocity and high inertia. This means that on long range shots 
the sabot may travel several miles down range and land with a high impact velocity. 
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Depending on the location of friendly forces and neutral populations along the line of 
fire, where the sabot conies down could be as important as where the projectile impacts. 


Front 


Figure (3) Nose and Body Configurations 
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Two methods for providing flat electrical contact areas on the projectile are 
considered. Figure (3) above illustrates both of these methods. The first method is to 
slice a cylindrical projectile lengthwise down the middle and insert a rectangular wedge 
of appropriate width with the front of the wedge matching the curvature of the nose. The 
second option for providing the electrical contact surface is to shave down opposite sides 
of a cylindrical projectile until a sufficiently wide flat area is planed out. This reduces 
cross sectional area, but lengthens the projectile for the same explosive payload. Neither 
has an intuitive advantage over the other so both are modeled in the configuration 
comparison program and a choice will be made by a trade study of the output. 
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The width of the electrical contact area must be considered. Making the electrical 
contact surface smaller makes the body more cylindrical,, which should reduce the 
separation and turbulence of the airflow around the body and therefore has aerodynamic 
benefit. A smaller contact area, however, means that the electric current has to go 
through a smaller area, which at the high currents required for rail gun operation leads to 
vaporization and spalling of the projectile's flat surface. Making the contact surface 
wider has the opposite effects of increasing separation and turbulence of airflow around 
the body and reducing current density. The two-inch wide contact area shown in Figure 
(3) is selected somewhat arbitrarily because it provides a sufficiently large electric 
contact area and the projectile remains fairly cylindrical. The width of the contact area is 
one of the physical parameters that can be changed by the user in the input section of the 
configuration comparison program to account for changing requirements. Since the 
contact surface extends the entire length of the body a narrow contact area may still 
provide a large electrical contact surface on both sides of the projectile. The drag 
equations in the configuration comparison program are written always assuming turbulent 
flow. Turbulent flow provides the worst drag which somewhat masks the effect of the 
flattened region. As specific parameters of the rail gun design are set and specific 
materials are chosen a smaller contact area may be possible. 

2. Nose Type 

There are several possible nose types for the projectile. The proposed maximum 
muzzle velocity of the rail gun design is Mach 8.9. The nose will experience extreme 
temperatures for the first several seconds after launch because of skin friction heating in 
the lower atmosphere and the shock front heating the air around the projectile. The high 
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temperatures will last until the projectile velocity is reduced, thereby reducing the 
strength of the shock, or the projectile reaches the upper atmosphere, where the lower air 
density will reduce friction. 

N. Kemp and F. Riddell [Ref. 3] show the heating of a projectile’s nose varies 
inversely with the square root of the radius of the nose’s point as seen in equation (2.1). 

Q = 20800(R n )" 1/2 (p/ p 0 ) V2 (u/ « c ) 325 ( 1 -h s k)h s i) (2.1) 

The other factors in the equation are environmental factors independent of the projectile. 
The temperature change of the nose is determined by equation (2.2). 

Mc p (dT/dt) = AQ (2.2) 

The factor AQ is the amount of heat being added by skin friction and the shock front 
minus the heat being radiated away by the nose to the atmosphere or conducted away 
from the nose to other parts of the projectile. A very pointed nose, meaning a small nose 
point radius, provides less drag, will experience more heating, and, because it has less 
mass over which to distribute that heat, risks melting and deforming. A very blunt nose, 
meaning a large nose point radius, has higher form drag, will experience less heat, and, 
because it has more mass over which to distribute that heat, is at reduced risk of melting 
and deformation. The rate at which the nose is heated and the thermal conductivity and 
emissivity of its constituent material determine the amount of mass required to absorb the 
heat without melting. If the nose is heated faster than it can radiate or diffuse the heat 
away if may melt or even vaporize. 

A tangent ogive with a fineness ratio, nose length divided by the maximum 

t 

diameter, of unity is used in the configuration comparison computer program. The 
fineness ratio is another of the physical parameters that the user can change in the input 
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section of the configuration comparison program. The user can select a fineness ratio to 
minimize the drag while maintaining heat capabilities for the conditions expected. A 
spherical nose, alternative to the tangent ogive, can also be modeled for comparison. 
Figure (3) on a previous page illustrates the two methods of providing the two-inch wide 
contact surface with both nose types and a two and one half -inch radius of body 
curvature. 

3. Component Material 

The material of which the projectile is made has important consequences on the 
eventual performance of the projectile. The material has to be a good conductor of both 
electricity and heat. A high-density material offers several advantages resulting in longer 
ranges and higher impact velocities. 

The biggest benefit of using high density material is the increased mass per unit 
volume. The drag force is determined by the geometry of the projectile and the 
environment in which the projectile is placed. The force on the projectile equals its mass 
times its acceleration. If the geometry of the projectile is not changed so the drag force is 
constant at a point in space, but the mass is increased then the acceleration, deceleration 
in the case of a drag force, must decrease as shown in Figure (4) on the next page. 
Therefore the effect of the drag force on the velocity of the projectile is inversely 
proportional to the mass of the projectile. The terminal velocity of the projectile is 
directly proportional to the mass of the projectile as calculated by equation (2.3). 

v t = (2mg/pS cs Ci>) 1/2 (2.3) 
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p is the ambient air density,Scs is the projectiles cross sectional area, and C D is the total 
drag coefficient. Another benefit of high-density material is the increased projectile 
penetration capabilities against hard targets. 


Figure (4) Acceleration vs Mess for 100 Newtons of Force 
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Tungsten is very dense and is already used for armor piercing projectiles. 
Tungsten is a good choice for this projectile and is used in the configuration comparison 
program for these reasons. The density of the component material can be selected by the 
user in the input section of the program. 

As with most thing there are always disadvantages to any option. For the sams 
reason the increased mass of a projectile made of tungsten resists the deceleration cause d 
by drag, it also resists the acceleration by the electric and magnetic fields in the rail gun, 
requiring higher currents and longer barrels. Shipboard portability is another major issue 
















for very heavy projectiles as the individual projectiles may no longer be handled by a 
single sailor. 

The disadvantages of large mass will probably drive an actual projectile to be 
made of a lighter alloy, possibly keeping an all tungsten nose for its penetration 
capabilities. There are limits, however, on how light a projectile can be made for a 
constant volume. As seen in Figure (4) the deceleration caused by drag will be greater 
for a lighter projectile and therefore the range and impact velocity will be little better than 
a conventionally fired projectile. 

The configuration comparison program can use the density of any material 
without effecting the accuracy of the modeling. The equations in the computer program 
that models the trajectory of the projectile, provided in Appendix (B) and fully developed 
in chapter four, are written assuming that the projectile remains supersonic all the way to 
impact. If the projectile becomes subsonic anywhere in the trajectory the calculations in 
the trajectory program diverge because several of the formulas are not defined for the 
subsonic regime. Therefore the projectile's component material must be sufficiently 
dense to keep the projectile supersonic until impact. 

4. Explosives 

A solid metal projectile at high velocities has a great deal of kinetic energy, but at 
long ranges where the projectile has slowed it may need to carry explosives to be 
effective. Against hardened targets, such as ships or concrete structures, the sudden 
deceleration of the projectile at impact transfers that kinetic energy into the target. 

Against soft targets, such as wooden buildings or troops, there is no mechanism to rapidly 
convert the energy of the projectile into destructive power at the target. The destructive 
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power is the impact energy divided by the time over which it is transferred. The shock 
wave of the projectile, if it is still at a high Mach number when it reaches the target, may 
cause some damage because of the pressure increase with the passing of the shock wave. 
The higher the Mach number of the projectile the stronger the shock will be. Figure (5) 
on the next page, computed with an equation presented by R. Zucker [Ref. 4] used to 
formulate the normal shock tables, shows the upper bound, since normal shocks are the 
strongest shocks, of the pressure ratio as the Mach number is increased. The pressure rise 
caused by the shock wave quickly dissipates with distance from the projectile and the 
obliquity of the shock front. The projectile itself will only punch a small hole through the 
walls of a wooden building or into the soft dirt like a present day armor piercing shell. 

Against either hard or soft targets at long ranges the projectile will decelerate to 
terminal velocity before impact and may not have much more kinetic energy than a 
convetional projectile. If the rail gun is to be a practical weapon it must be able to 
perform at least as well as the current weapons, such as the five inch gun it may someday 
replace. This means the proposed projectile will have to carry an explosive payload 
similar to a convetional projectile. Twelve pounds have been selected as the payload 
mass for the model because that is the mass of explosive material in a five inch shell. 

The explosive compound C4 is used in the computer model because the density of C4 is 
about average for modem high energy explosives. The explosive payload and the density 
of the selected explosive are both variables which may be changed by the user in the 
input section of the computer model. The length of the projectile is calculated from the 
volume required to contain the selected mass. Substituting a different explosive 
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Figure (5) Pressure Ratio for Passage of Normal Shock 



compound of similar density at the same explosive payload should only make a small 
difference in the length of the projectile. 

5. Casing Thickness 

The explosive for the projectile must obviously be encased within the projectile so 
the question of casing thickness arises. A casing thicker than a convetional projectile 
would provide more material to help absorb the expected aerodynamic heat transfer. The 
projectile's cross sectional area and surface area, however, will be increased for a 
constant explosive payload which increases the drag. The case thickness is yet another 
variable that a user can select in the model's input section. A half inch casing thickness 
was used in the model presented in the final section of this chapter. 







Additional options regarding materials emerge at this point. Using a single 
material for the entire projectile eases manufacture since the nose and body can be cast as 
a single unit, so the model uses a single material, but the computer model can easily be 
altered to make the nose, end plate and the casing out of different materials. This is 
significant because the nose is not in contact with the conducting rails rails so its electric 
conductivity isn't very important and the casing will be the main conductor of electricity 
in the barrel so its conductivity is of great importance. 

B. DIMENSIONS AND FORCES 

This section details the development of the equations for the configuration 
comparison program that define the physical dimensions of the projectile's nose and 
body. Many of the calculations are complicated by keeping the physical dimensions as 
variables to allow the flexibility in the computer model discussed in the previous section. 

1. Nose 

The mass and surface area of the nose are the important parameters for 
determining the the effects of the aerodynamic properties in the equations of motion. 
Refering back to Figure (3) will help the reader visualize the nose configurations. To 
determine the mass, the volume of the nose and the mass density of its component 
mate rial are required. Once the volume is determined the mass can be calculated by 
multiplying the volume and the density. 

The geometry of the wedge insert configuration is simpler, making the required 
parameters easier to calculate. The volume of a half sphere and conical ogive are known 
formulas, the formula for the ogive can be found in Appendix A of Reference (1). The 
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inserted wedge for the spherical nose is just a half cylinder of the same radius and desired 
width. Using r as the radius of the spherical nose and h as the width of the contact area in 
equation (2.4). 

NVs/ w = (2/3)ur 3 + (l/2)7tr 2 h (2.4) 

The calculation for the ogive wedge configuration is a little more involved. Calculus 
provides a convenient solution for the area of a section of a circle, shown in Figure (6) 
below. To get the volume we multiply the planform area that matches the nose by two 
and the desired width. Equation (2.5) uses L as the nose length, R as the radius of 
curvature of the nose, and r and h as above. 

NV o/w = 7t[L(2R 2 -2Rr+r 2 ) - L 3 /3 - (R-r)L(R 2 -L 2 ) 1/2 - (R-r)R 2 sm 1 (UR)] + 

[R 2 (ti/ 2-0) - (R-r) 2 cot0]h - (2.5) 

Figure (6) Calculating Radius of Curvature and Planform Area for Ogive 

r = Radius of Projectile 
L = Nose Length 



solving the integrals 


|[R 2 (^-0)-(R-r) 2 cotG] 

The volume of the shaved side projectile configuration is easily understood in 
concept, but mathematicaly complicated to compute. The concept is simply to determine 
the volume for the complete nose shape and subtract the volume of the piece that is 
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shaved away. The removed volume is determined by solving the triple integral in 
equation (2.6). 


h/2 n/2 R 

RV = J / J r'dr’de’dz (2.6) 

0 cos'^fz/R) (R 2 -fV) 1/2 csce’ 

The term f is the nose fineness ratio. Completing the integration with variables for the 

limits of integration requires an approximation. Simpson's Rule is used to get the 

solution of the integral above for equation (2.7): 

RV = R 2 hit/8 - (R 2 h/4)cos' 1 (fh/2R) - (R 3 /2f)(l-(fh/2R) 2 ) 1/2 + R 3 /2f - (2.7) 

(Rfh/48) [(h/2)csc(cos _1 (hf/ 8R)) + (h/2)csc(cos' 1 (hf , 4R)) + 

(3h/2)csc(cos' 1 (3hf8R)) + (h/2)csc(cos' 1 (hf/2R))] + 
(f 3 h/48R) [(h 2 /12 8)csc(co s' 1 (hf78R)) + (h 2 /32)csc(cos' 1 (h£ , 4R)) + 

(27h 2 /128)csc(cos‘ 1 (3hf / 8R)) + (h 2 /8)csc(cos' 1 (hf/2R)) 


Multiply ing this volume by four gives the total volume lost from the nose by planing the 
flat sides on to the projectile. A spherical nose is the limiting case of the ogive nose 
where the fineness ratio equals one half, the radius of curvature of the nose and the length 
of the nose equal the radius of the projectile. The volume equation for the shaved side 


spherical nose configuration is equation (2.8). 

NVs/s = (2/3):ir 3 - R 2 lm/2 + R 2 hcos' 1 (h/4R) + 4R 3 (l-(h/4R) 2 ) 1/2 - R 3 + (2.8) 

(Rh/24)[(h/2)csc(cos' 1 (h/16R)) + (h/2)csc(cos' 1 (h/8R)) + 

(3h/2)csc(cos' 1 (3h/16R)) + (h/2)csc(cos' 1 (h/4R))] - 
(h/96R)[(h 2 /128)csc(cos' 1 (h/16R)) + (h 2 /32)csc(cos' 1 (h/8R)) + 

(27h 2 /128)csc(cos‘ 1 (3h/16R)) + (h 2 /8)csc(cos’'(h/4R)) 

The volume for the shaved side ogive nose configuration is shown in equation (2.9): 

NVo/s = 7t[L(2R 2 -2Rr+T 2 ) - L 3 /3 - (R-r)L(R 2 -L 2 ) 1/2 - (R-r)sin'‘(L/R)] - (2.9) 

R 2 hrc/2 + R 2 hcos' 1 (fh/2R) + (2R 3 /f)(l-(fh/2R) 2 ) 1/2 - 2R 3 /f+ 

(Rfh/12) [(h/2)csc(cos‘ 1 (hf/8R)) + (h/2)csc(cos‘ 1 (hf/4R)) + 
(3h/2)csc(cos' 1 (3hfr8R)) + (h/2)csc(cos‘ 1 (hfr2R))] - 
(f 3 h/12R)[(h 2 /128)csc(cos' 1 (h£ / 8R)) + (h 2 /32)csc(cos 1 (hfr4R)) + 
(27h 2 /128)csc(cos' 1 (3hfr8R)) + (tf^csctcos^Oif^R)) 
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The surface area of the various nose configurations, needed to determine skin 
friction drag. Once again the wedge insert configuration surface areas were easier to 
calculate. The spherical nose with wedge insert is the surface area of a half sphere plus 
the surface area of half a cylinder of matching radius and desired electrical contact area 
width. The formula for the spherical nose with wedge insert is given in equation (2.10): 

Sns/w = 27tr 2 + 7trh (2.10) 

The equation for the surface area of the ogive nose is a surface of revolution and is 
determined by solving the integral in equation (2.11): 

2fr 

So = 1 27i[(RV) 1/2 - (R-r)][l+(y 2 /(RV)] 1/2 = 2;r[2frR - (2.11) 

0 R(R-r)sin' 1 (2fr/R)] 

The surface area of the inserted wedge is the desired electrical contact area width and the 

arc length, which is determined from the radius of curvature of the nose and 0 as 

calculated in Figure (6), doubled to account for the upper and lower surface of the wedge. 

The surface area of the ogive nose with wedge insert is the surface area of the ogive nose 

plus the surface area of the wedge insert, shown in equation (2.12). 

S„ 0 /w = 2u[2frR - R(R-r)sin’ 1 (2fr/R)] + 2hR0 - (2.12) 

The surface area of the shaved side configuration is calculated in much the same 
manner as the volume was determined. The surface area of the removed section is 
subtracted form the surface area of the whole surface, but now the surface area of the 
remaining flat surface must be added back to the surface area. For the shaved side 
spherical nose one can use the symmetry as we did in the volume calculation. The 
removed surface area can be considered a surface of revolution and determined by 
solving the integral in equation (2.13). 
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.5h 

Jo 27ty(l+(dy/dx) 2 ) 1/2 dx = 27t[{-.5r ln(.5h + (rf.25h 2 ) I/2 )} + 

{.25h(r+.25h 2 ) 1/2 } + {.5r ln(r)}] (2.13) 

The integrand in equation (2.13) is y = (r^x 2 ) 172 and the limits of integration are from the 
centerline of the projectile to the edge of the flat electrical contact surface. Each of the 
resulting flat areas left on the nose is a half circle with a radius of one half the electrical 
contact area width. Collecting all the terms for the surface of the half sphere, the surface 
of the removed portion, and the surface of the remaining flat spots results in equation 
(2.14). 

- 2%? - 27i[{-.5r ln(.5h+(r+.25h 2 ) 1/2 )} + (2.14) 

.25h(r+.25h 2 ) 1/2 } + {.5r ln(r)}]+7i(h/2) 2 

The shaved side ogive nose is harder to calculate because of the lack of symmetry. The 
removed surface area for the shaved side ogive nose is four times the value of the double 
integral solved in equation (2.15). 

h/2 (tc/ 2 - cos'^fe/r)) 

J | R d0 dz = R[7ih/4 - (h/2)cos 1 (fh/2R) - (2.15) 

0 0 (R/f)(l-(fh/2R) 2 ) 1/2 + R/fJ 

The flat area left on the nose when the sides are shaved down is once again shaped like 

the area in Figure (6). The surface area of the ogive nose with shaved sides configuration 

is the combination of the surface area of the ogive nose minus the surface area of the 

removed section plus the surface area of the flat spots left on the nose when the shaved 

area is remove. The ogive nose with shaved sides surface area formula is given by 

equation (2.16). 
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(2.16) 


Sws = 2?t[2frR - R(R-r)sm 1 (2fr/R>] - K[nh- (2h)cos 1 (fh/2R) - 
(4R/f)( 1 -(fh/2R) 2 ) 1/2 + 4R/f] + 2[R 2 (7t/2-cos' 1 (h/R» - 

(R-h/2) 2 cot(cos" 1 (h/R))] 

2. Body 

The cross sectional area, length, mass, and surface area are the important physical 
characteristics of the projectile body. The cross sectional area of the wedge insert 
configuration is the area of a circle who's radius is the radius of the projectile plus the 
area of a rectangle who's height is twice the radius and who's width is the width of the 
contact area, as shown in equation (2.17). 

Scsw = Ttr 2 + 2rh (2.17) 

The cross sectional area of the shaved side configuration is the area of a circle minus four 
of the planform areas that were calculated for the shaved side spherical nose, shown in 
equation (2.18). 

S oss = Ttr 2 - 2[r 2 (7i/2-0) - (r-rfcotO] (2.18) 

The length of the projectile's body depends on how much explosive is required. 
First the cross sectional area of the explosive must be determined for the desired 
projectile radius and skin thickness. The same equations as were used for the body's 
cross sectional area can be used to determine the explosive's cross sectional area by 
reducing the radius used in the equation by the desired skin thickness. The length can 
now be determined by dividing the desired explosive payload mass by explosive cross 
sectional area and the density of the chosen explosive. Equation (2.19) uses me as the 
mass of the explosive, p e as the density of the explosive, and S ecs as the appropriate 
explosive cross sectional area for the wedge or shaved configuration. 

Lb Hle/peSecs (2.19) 


22 


The mass of the projectiles body is the sum of the explosive mass plus the mass of 
the casing mat erial. To determine the mass of the casing it is easiest to find the volume 
of the material by subtracting the explosive cross sectional area from the projectile's cross 
sectional area and multiplying by the length of the explosive and density of casing 
mater ial Equation (2.20) uses the variables from above and Spcs for projectile cross 
sectional area, p c for the density of the casing material. Lb for the body length, and mb for 
mass of the body. 

nib = ®e LbPc(Spcs ” Secs) (2.20) 

Of course the length and cross sectional areas must all be for the same wedge or shaved 
configuration. 

The body's surface area for the wedge insert configuration is, again, easily 
determined. It is the surface area of a cylinder plus the upper and lower contact surface 
which is the width of the contact area times the length of the body as shown in equation 
( 2 . 21 ). 

S bw = 27trLbw + 2hLbw (2-21) 

The shaved side configuration requires reference back to Figure (3). In order to get the 
arc length to determine the surface area the figure should be viewed on its' side. For this 
case we know that RsinG equals half the contact area width where R is the radius of the 
projectile. Solving for 0 and multiplying by the radius determines half the arc length for 
one side of the projectile. Multiplying the arc length by four and by length of the shaved 
side body gives the surface area of the curved portions of the projectile body. The flat 
surfaces on the projectile are once again the contact area width times the body length for 
the shaved side configuration. Combining the surface areas for the shaved side 
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configuration we get equation (2.22), the total surface area for the body in this 
configuration: 

S bs = 4rL bs sm'(h/r) + 2hL bs (2.22) 

3. End Plate 

The end plate is assumed to be the same thickness as the casing, although this can 
be changed in the programs if desired. Since the thickness is set the only parameters to 
find are the surface area and mass. The surface area can be easily calculated by taking 
the equation for the body surface area for the matching configuration from the previous 
section and replacing the body length with the skin thickness. The mass is just as easily 
determined by multiplying the cross sectional area of the matching body configuration by 
the skin thickness and component material density. 

4. Drag 

The form and. skin friction drag for the different projectile configurations are 
computed using Reference (1) to compare the drag characteristics of each. The dynamic 
head is required for drag calculations and is computed in the model using standard sea 
level air density and the selected muzzle velocity with equation (2.23): 

q = (l/2)pv 2 (2.23) 

The skin friction drag is computed using a generic high velocity skin friction drag 
coefficient form Reference (1) of .002. The generic skin friction drag coefficient is an 
acceptable approximation because the large velocity of the projectile is the Ramp for all 
cases and dominates the numerator of the Reynolds number calculation, shown in 
equation (2.24) used to determine an exact skin friction drag coefficient, shown in 
equation (2.25). 
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Re = pvL il\x (2.24) 

C f = .455 [Logio(Re)]' 2 ' 58 (2.25) 

The skin friction drag, equation (2.26), is the skin friction drag coefficient multiplied by 


the total surface area and the dynamic head. 

C Df =C f Sq (2.26) 

The form drag calculation has a few more intermediate steps. The form drag coefficient 
multiplied by the projectile's cross sectional area, Spcs, and the dynamic head previously 
calculated determines the form drag in equation (2.27): 

Cdf = C F (Spc S )q (2.27) 

The form drag coefficient is calculated with equation (2.28): 

C F = (Ap/q) {1 -[(196(L n /d) 2 -16)/( 14(M+18)(L„/d) 2 )]} (2.28) 

The conical shock pressure coefficient is calculated using the empirical formula 
given in equation (2.29). 

Ap/q = [.083+(.096/M 2 )](<t/10) 1 ' 69 (2.29) 

The semivertex angle, a, of a reference cone corresponding to the ogive nose is 
required to determine the pressure rise across a conical shock and is detennined in 
equation (2.30) by t aking the arctangent of one divided by twice the fineness ratio. 

<y = tan _1 (l/2f) (2.30) 

B. RESULTS 

The configuration comparison computer program computes the physical data for 
the projectile over a range of projectile radii selected by the user for the different 
projectile configurations and plots the results for comparison. Figures (7) through (12). 
on the next several pages presents these plots for one computer program run with the 
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Form Drag/io 5 (kg-m/a 2 ) 


Figure (7) Form Drag vs Skin Friction Drag at Sea level, V0=3000 km/s 
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Figure (12) Skin Friction Drag vs Width at Sea Level, V0=3000 km/s 
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projectile radius ranging from two to six inches. The majority of the inputs and the plots 
are in English units, but the model can display the plots in SI units by minor 
modifications to the variables in the plot statements. 

If the only consideration is to maximize the range and impact energy of the 
projectile, then minimizing the total drag will optimize the projectile's performance. The 
form drag is two orders of magnitude larger than the skin friction drag in Figure (7) so 
the minimum total drag is along the bottom of the plot. 

Examining Figure (8) shows that the minimum form drag is the shaved side ogive 
configuration with a four-inch diameter. Refering to Figure (9) the four inch wide shaved 
side ogive projectile is forty one inches long. That is not a very practical length. Figure 
(10) or Figure (11) shows that the projectile would weigh one hundred and seventy-eight 
pounds. 

The five inch diameter shaved side ogive projectile has only slightly higher form 
drag as shown in Figure (8) and substantially lower skin friction drag shown in Figure 
(12). The length of the five inch shaved side ogive projectile in Figure (11) is twenty 
four inches. Twenty four inches is a far more practical length for many reasons. A 
twenty four inch projectile is easier to store, will have less longitudinal stresses in flight, 
and the explosive material is easier to detonate when it is packed together rather than 
stretched out. Refering to Figure (10) or Figure (11) the five inch shaved side projectile 
would weigh one hundred and forty-four pounds. 

Using the same figures as the previous two cases, the six inch diameter shaved 
side ogive projectile is shorter still at only fifteen inches in length, but the form drag is 
beginning to get appreciably larger. The skin friction drag for the six inch shave side 
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ogive projectile is actually the minimum, but the form drag increase is substantially larger 
than the skin friction drag decrease. The six inch shaved wide projectile would weigh 
one hundred and forty-seven pounds. The increased drag without increased weight 
would mean a decrease in performance. 

Many of the assumptions made in defining the inputs to run the model were 
selected with the idea of replacing a conventional five inch gun with a rail gun, so it 
should be no surprise that a modified five inch projectile is the apparent best choice. As 
other requirements are specified the plots can be used to make similar comparisons 
between the parameters to best meet the desired performance. There are additional plots 
produced by the configuration comparison program that weren't needed for the above 
analysis, but may be useful for other comparisons. Many more parameters are calculated 
in the model than are plotted. A user can change or add plot commands to get the plots 
necessary to make the comparisons they desire for the requirements they are given. 
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III. STABILITY 


Projectiles fired over any appreciable distance require some form of stability to 
keep from being tumbled by minor variations in wind speed i.e., wind shear. Projectiles 
fired from conventional guns are spun in the barrel so they have angular momentum to 
provide stability. The proposed rail gun barrel, being rectangular, will not impart this 
an gular momentum to its projectile. Fins are the most practical means of providing a 
stabilizing force. 

A stability program, also written in MATLAB and provided in Appendix (C), 
models the stability of a projectile with fins. The shaved side projectile configuration 
appears better able to provide long range and high impact velocity than the wedge insert 
configuration because of the shaved side's lower drag. Therefore, the shaved side 
configuration is modeled in the stability program. The equations for the ogive nose are 
used in the model, but if a spherical nose is desired it is simply an ogive with a fineness 
ratio of one half. 

A. FINS 

A variety of ideas were initially considered for attaching the stabilizing fins to the 
projectile. Most of the ideas considered for attaching fins were concerned with whether 
or not there was sufficient room in the barrel to accommodate fins large enough to 
provide stability. Attaching large fins requires that they be at least partially retracted 
while in the barrel of the rail gun and that they deploy after the projectile has left the 
barrel. Folding fins add complexity and, as such, cost to the over all system. 

Selecting short span fins with long chord, as shown in Figure (13) on the next 
page, it will be shown will provide the stability without adding the complexity and cost of 
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Figure (13) Short Span Long Chord Fin 
Planform View 



Front View 

W- b-H 


a folding fin configuration. A cruciform configuration with the fins extending to the 
comers of the rectangular barrel is probably the simplest configuration and is what was 
shown in Figure (1). Symmetric biconvex fins are used in the stability program because, 
as per Reference (1), they are easy to manufacture and offer the best structural strength. 

The dimensions for the nose and body are calculated with the shaved side ogive 
nose equations developed in the previous chapter. The dimensions of the fins are 
calculated by the model using several inputs from the user and values from the nose and 
body calculations. 

1. Span and Chord 

The span of the individual fins is calculated in equation (3.1) by subtracting the 
radius of the projectile from the hypotenuse of a triangle whose sides are half the 
projectile diameter at the shaved surfaces. 

b = [(r-r’) 2 +(r-r') 2 ] 1/2 - r (3.1) 
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Calculating the span in this way keeps the fins orthogonal to each other and keeps their 
length correct to just touch the barrel wall regardless of selected projectile radius. Many 
of the aerody nami c calculations assume the fins in pairs without the body present to form 
a clipped delta whose span is just twice the individual fin's span. 

The root chord is determined as a user selected percentage of the total body length 
calculated for the selected body radius. The tip chord is calculated in equation (3.2) by 
multiplying the span by the tangent of the sweep angle of the fin's leading edge. A, and 
subtracting the result from the root chord of the fin. 

C t = Q - btanA (3.2) 

The selected chord length has two effects in the calculations of the other physical 
dimension. The longer the chord length, the shorter the span can be while still providing 
stabilization of the projectile. The thickness of the fin is determined by multiplying the 
chord by a user inputted thickness to chord ratio, t. The longer the chord length, the 
thicker the fin becomes which increases form drag. Users are cautioned to remain 
conscious of the thickness of the fins as form drag caused by unnecessarily thick fins may 
severely degrade the performance of the projectile. 

The sweep angle is set in the stability program to keep the leading edge of the fin 
forward of the Mach cones created at the point of the delta, which keeps the leading edge 
supersonic until the freestream becomes transsonic. If a user is sure that the projectile 
will never slow to transsonic velocities, a larger sweep angle may be used without the 
calculations, such as for the fin form drag coefficient, diverging. The sweep angle is set 
in the fins section of the program rather than in the input section to prevent inexperienced 
users from getting in trouble by making a poor selection of sweep angle. 
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2. Surface Areas 


The planform area of a fin pair is used as a reference area in many intermediate 
calculations relating to the fins, all the final quantities are referenced to the maximum 
body cross sectional area. The simplest way to calculate the planform area is to break up 
the fin pair into a pair of triangles and a rectangle. The area of the two triangles is the the 
root chord minus the tip chord multiplied by the span. The rectangle is the tip chord 
multiplied by the span of a fin pair. Equation (3.3) combines the terms. 

Sf=(C r -C t )b + 2C t b (3.3) 

Since the fin configuration is a clipped delta wing some of the calculations are 
based on what percentage of the full delta wing is present, so the planform of the full 
delta wing is required. The planform of the full delta wing is the root chord multiplied by 
half the span, as shown in equation (3.4) 

S d = C r b (3.4) 

3. Volume and Mass 

The thickness of the fin is determined from a thickness to chord ratio, t, so as the 
chord lessens from root to tip so does the fin thickness. This makes determining the 
volume of the fins very difficult. Since the fins do not have a large sweep angle and do 
have a very short span using the surface area of a rectangular fin with the mean 
aerodynamic chord of the actual fin offers a very close approximation to the surface area 
of the actual fin. The mean aerodynamic chord for a clipped delta wing is calculated with 
equation (3.5): 

C™ = (2/3)0-{ [ 1 +(Ct/Cr)+(Ct/C r ) 2 ]/[ 1 +(Ct/Cr)]} (3.5) 
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The volume can now be approximated by determining a rectangular fin cross 
sectional area at the mean aerodynamic chord and multiplying it by the span. To 
determine the fins cross sectional area the equations from Figure (6) can be used once 
again. The symmetry of a rectangular biconvex wing allows the calculation for a single 
section to be multiplied by four to get the total volume. The volume of the fins and the 
intermediate steps required to calculate it are shown in equations (3.6), (3.7), and (3.8). 


Rf- [(-5C m a) + (.5TC ma ) ]/x C ma 

(3.6) 

Of = COS'^.SCma/ Rf) 

(3.7) 

Vf = 4b[(Rf) 2 (7t/2 - Of) - (Rf - .5TC ma )cot0f] 

(3.8) 


The mass of the fins can then be calculated in equation (3.9) by multiplying the volume 
bytwo, for the two fin pairs, and the density of the component material. 

m f = 2 Vfpc (3.9) 

B. MATHEMATICAL MODEL 

The physical dimensions calculated are now use to determine aerodynamic 
coefficients. The combination of the aerodynamic coefficients will be used to determine 
if the projectile is stable. 

1. Center of Gravity 

The center of gravity is required to calculate the stability. Since the projectile is 
symmetric along its length the simplest method of calculating the center of gravity is to 
add the first moments of all the components, which will all be on the center axis, and 
divide by the total mass. Placing the origin at the point of the nose the first moment 
equations for the nose, body, end plate, and fins become equations (3.10), (3.11), (3.12), 
and (3.13) respectively. 
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I„ = .707(2fr)m„ 

(3.10) 

lb — [. 5L b 4*2fr] [p c L b (S pcs - S ecs )+me] 

(3.11) 

Iep = [. 5(st)+Lb+2fr] [pcCstXSpcs)] 

(3.12) 

If = [-5 C ma +( 1 -(%+. 1 ))L X ] [2nif] 

(3.13) 


The percent sign is the root chord percentage of body length and Lx is the total projectile 
length. 

2. Graphical Approximations and Miscellaneous Factors 

Many aerodynamic coefficients and factors are determined empirically and given 
in graphical form. The four interference factors from Reference (2) and the nose/body 
lift coefficient from Reference (5) are presented graphically and are required for several 
of the calculations. The user of the stability program would have to calculate many of the 
projectile's physical parameters to read the graphs and find the correct inputs to the 
program everytime a parameter, such as projectile radius or root chord percentage, was 
changed. Having the user do many of the calculations defeats the purpose of having a 
computer model. Computers, unfortunately, cannot read graphs so it became necessary 
to determine equations that the computer could use to calculate the values of coefficients 
and factors. Appendix (D) has a series of short MATLAB programs written to determine 
curves that match the graphs from the above references. The plots from these MATLAB 
programs are provided in Appendix (D) with several points, corresponding to actual 
points from the references, marked for accuracy comparison. Reviewing these plots will 
reveal that the approximations are seldom more than a couple of percent off. 
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There are several calculated aerodynamic quantities that are used in the equations. 
The aspect ratio calculated in equation (3.14) is the square of the span divided by the 
planform area of the fin. 

AR = b 2 /S f (3.14) 

The Mach angle p is the angle of the semivertex of the Mach cone and is determined in 
equation (3.15) by taking the inverse sine of the inverse of the Mach number. 

p = sin' 1 (1/M) (3.15) 

The final quantity is p, which is the square root of the Mach number squared minus one, 
shown in equation (3.16). 

p = (M 2 - 1) 1/2 _ (3.16) 

3. Center of Pressure 

The center of pressure of the nose/body is also a graph determined quantity from 
Reference (5), but the plot is a fixed straight line so long as the body length is at least 
twice the nose length. Similar to the nose/body lift coefficient, shown in Appendix (4), 
the center of pressure graph is divided into two separate lines at the point where p divided 
by the nose fine ness ratio equals one. An if/then statement in the stability program 
determines when P divided by the nose fineness ratio is less than or equal to one and then 
calculates the two quantities with equations (3.17) and (3.18): 

CNan - -.28[(p/f)-1.9] 2 + 3.58 (3.17) 

XcPn = CG - [(.125L T p/f)+.15] (3.18) 

The center of pressure is referenced to the center of gravity. If P divided by the nose 
fineness ratio is greater than one the program uses equations (3.19) and (3.20): 

CN an = -3.4[(frP)-.9] 2 + 3.4 (3.19) 
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Figure(14) Mach Cones on Fin Planform 


Crossed 



Hie effected area is 
represented by the 
striping. 


Uncrossed 



Xcrn = CG - [(-. 125L T f/P)+.425] - (3.20) 

The area of the fin effected by the Mach cones created at the point of the delta and 
at the fin tips is defined by the p calculated above. The effected surface area will 
increase as the Mach number decreases and decrease as the Mach number increases. The 
effected area has a lift coefficient that is half the coefficient for the fin surface in the 
uneffected area and is thus an important intermediate step in calculating the total lift 
coefficient and center of pressure location of the fins. Calculating the effected surface 
area is also done by dividing the calculation into two regions. If the span is short and the 
root chord is long, as in Figure (1), the Mach cones will cross somewhere on the fin's 
surface, but if the root chord gets shorter and the span longer or the Mach number gets 
very large the Mach cones may not cross until after the fin's trailing edge, as shown in 
Figure (14) above. In order not to double count the effected surface area another iftthen 
statement is used. Whether the Mach cones have crossed before the trailing edge is 
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determined by adding the base at the fin's trailing edge of the triangles formed on the 
planform surface by the Mach cones and comparing to see if that exceeds the span. If the 
Mach cones have crossed the unaffected planform area is computed and subtracted from 
the total fin planform area with equation (3.21). 

S e ff = Sf- {[(1.1452tan(7i/2-A-p))/(tan(7i/2-A-ji)+tan(A-Tr/2-p))] 2 + (3.21) 

[(1.1452tan(A-7r/2-p))/(tan(n/2-A-p)+tan(A-7i/2-p))] 2 }tan(7i/2-A-p) 

If the Mach cones do not cross the calculation of the effected planform surface area is 

simp ly the addition of the surface area of the two triangles formed by the Mach cones 

shown in equation (3.22). 

S e ff = Cr 2 tanp + C 2 tanp (3.22) 

The normal force coefficient of the fins can now be determined. If there were no 
Mach cones the normal force coefficient would simply be four divided by the quantity p. 
The effected area in the Mach cones, however, only has a normal force coefficient that is 
half the value for the uneffected area, two divided by the quantity p. To account for the 
reduced lift coefficient in the effected area half the lift coefficient multiplied by the 
percentage of total fin planform area effected is subtracted from the total as shown in 
equation (2.23). 

CN af = (4/p) - (2Seff/pS f ) (2.23) 

The center of pressure for a full delta wing is simply two thirds the root chord. 
Since the fins being used are basically clipped delta wings the center of pressure is 
moved forward and related to what percentage of the full delta wings area is actually 
present. The center of pressure of the fin is also effected by the Mach cones. The Mach 
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cone effected area only counts as half its actual area as a percentage of the total delta 
wing area. The center of pressure for the fin is calculated by equation (3.24). 

Xcrf- (2/3)C r (Sf/S de „a - .5Seff/S d e,ta) (3.24) 

The center of pressure above is referenced to the leading edge at the root. In order to 
combine the fin center of pressure with the nose/body center of pressure it must be 
referenced to the same point. The nose/body center of pressure was easily referenced to 
the center of gravity since the tip of the nose is always at the origin of the coordinate 
system. The position of the fin’s leading edge, however, will change everytime the root 
chord percentage or the radius of the projectile is changed. Foreshadowing the drag 
calculations of the next chapter, there will be a boattail on the projectile and the trailing 
edge of the fin will be at the start of the boattail. The boattail is set at ten percent of the 
body length-so the leading edge of the fins referenced to the nose of the projectile will be 
at the total length of the projectile minus the root chord percentage plus ten percent. 
Subtracting the center of gravity location from the fin leading edge location will move the 
reference to the center of gravity frame. The center of pressure calculation for the fins 
becomes equation (3.25). 

X CP f = (2/3)C r (Sf/SdeIta - .5S eff /Sdeita) - [(1-(%+.1))L t -CG] (3.25) 

This center of pressure that relates the interference between the nose/body and 
fins and is found in graph (15) of Reference (2). As long as the tip chord does not 
become substantially shorter than the root chord and the span does not become 
substantially larger than the diameter of the projectile body equation (3.26) using a linear 
approximation of the graph will work. 

X C Pb = -[ Cr(.9988p(AR)+.35)] - [(1 -(%+.!))L T -CG] (3.26) 
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Graph (15) of Reference (2) once again is referenced to the root chord leading edge so the 
second factor is used to rereference the position to die center of gravity. If a 
configuration outside the limits above is desired the user will have to recalculate the 
above approximation using Reference (2) to meet their requirements. 

4. Pitching Moment Coefficient 

All the calculations above have been leading to the most important calculation. If 
the pitching moment coefficient, Cm a , is negative there will be a restoring force acting 
against angle of attack changes to keep the projectile stable. The slope of the pitching 
moment coefficient in equation (3.27) is the combination of the quantities, most 
imp ortantly the center of gravity referenced center of pressures, calculated above. 

CMa = {CN a Xcp„ + [CNaf (X CPA, K B w + x cpf Kws)S/2r]}/2r (3.27) 

C. RESULTS 

The stability program, much like the configuration comparison program from the 
previous chapter, plots many different quantities against Mach number. Many more 
combinations are possible, should the user desire. The important plot for this thesis is the 
plot of the pitching moment coefficient as the Mach number is varied. 

The projectile from Figure (1) with a long root chord, fourteen inches, and short 
span, about three quarters of an inch, and a five inch shaved side body with an ogive nose 
was modeled. T his model produced the results presented in Figure (15) on the next page. 
Figure (15) shows that the slope of the pitching moment coefficient is indeed negative 
through the Mach region of interest, so there will be a restoring force against any angle of 
attack changes. An important distinction should be made here. The fins will help 
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Figure (15) CM alpha 



prevent the projectile from being tumbled by winds encountered in flight, but they will 
not prevent the projectile from being offset laterally by the wind, just like with 
conventional projectiles. The next chapter examines the trajectory of the projectile and 
depending on the angle of fire it is possible that the projectile may pass through the high 
winds of the jetstream. The speed of the jetstream winds, normally no higher than one 
hundred and fifty meters per second according to Reference (6), is small compared to the 
velocity of the projectile, normally no lower than one thousand meters per second at 
altitude. If the projectile were moving at its slowest and at ninety degrees to the 
maximum jet stream wind the angle of attack would only be eight degrees. The angle of 
attack change will be resisted by the fins, but the impact point of the projectile may be 
moved substantially. 
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rv. PERFORMANCE 

The previous chapters have provided the physical dimensions of a stable 
projectile that meets the requirements of the rail gun concept. This chapter reviews the 
development of a program, written in the C computer language and provided in Appendix 
(B), to model the performance of the projectile. 

A. PROJECTILE 

Most of the dimensions of the projectile and fins are calculated with the equations 
from the configuration comparison and stability programs; translated into C for the 
performance program. There are, however, a couple of quantities not previously needed. 
The wetted area of the fins is used for the skin friction drag calculation and the base area 
is used for the base drag calculation. 

The wetted area of the fins is determined by using the same assumptions that were 
used in determining the volume of the fins. The arc length at the mean aerodynamic 
chord is calculated using the radius of curvature of the fins and the angle as determined in 
Figure (6). The wetted surface area of the fins is the product of the radius of curvature, 
the angle subtended, the span, and by four to account for the four such surfaces that make 
up each fin. Equation (4.1) shows the formula using the same notation as in the fin 
volume calculation of the previous chapter. 

Swf-4bRf0 f (4.1) 

The base area of the projectile is the area at the very end of the projectile. The 
base area is reduced from the projectile’s maximum cross sectional area because of the 
boattail. The boattail is set at eight degrees and as ten percent of the total body length. 
Reference (1) shows that these parameters provide the best performance at lower Mach 
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numbers which is where the projectile will spend the majority of its flight. The radius of 
curvature and the height of the rectangular area will be reduced by one tenth the total 
length multiplied by the tangent of eight degrees as shown in equations (4.2), (4.3), and 
(4-4). 


r ba = r-.014054L T 

(4.2) 

0 ba = 7i/2 - sin' I (h/2r) 

(4.3) 

H ba = r ba cos(sin' 1 (h/2r)) 

(4-4) 


Using equation (2.18) used for the projectile cross section, but with the reduced quantities 
from above, the surface area of the base can now be determined with equation (4.5). 

S ba = 2[(r ba ) 2 sin *(h/2r) - (r ba ) 2 cos 2 e ba cot0 ba ] - 

2hr ba cos(sin' 1 (h/2r)) . (4.5) 

B. ATMOSPHERE 

The drag force on the projectile is dependent on the projectile’s environment. As 
the projectile moves through the atmosphere that environment changes, so the 
atmosphere, i.e., the environment, must be modeled. A standard atmosphere can be 
modeled as a series of linear temperature gradients and isotherms with increasing altitude 
as shown by reference (7). The temperature gradients, in degrees Kelvin are all of the 
form given in equation (4.6). 

T = T 0 + ay/1000 ’ (4.6) 

Where To is the temperature at the bottom of the gradient in degrees Kelvin, a is the slope 
of the gradient in degrees Kelvin per meter, and y is the altitude in meters. As an 
example the To is two hundred and eighty eight degrees Kelvin and the slope is minus six 
and a half from zero to eleven kilometers. The temperature is needed to determine the 
speed of sound and air density at each altitude. 
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The speed of sound, equation (4.7), at any altitude is the square root of gamma 
multiplied by the gas constant, R, for air and the temperature: 

a = (yRT) m (4.7) 

The Mach n umb er, used in the drag calculations below one hundred and five kilometers, 
is determined by dividing the projectile's velocity by the local (at altitude) speed of 
sound. 

The skin friction drag of the projectile at altitude is directly proportional to air 
density at that particular altitude. The air density for any altitude is calculated in equation 
(4.8) as a percentage of sea level air density with the following equation: 

p a = pao(T/288)- |yg/R(T ' 288)] ' 1 (4.8) 

C. DRAG 

The performance program uses six drag coefficients in determining the trajectory 
and impact velocity of the projectile. The skin friction drag and form drag coefficients 
for the nose/body and fins, and the base drag and boattail drag coefficients are calculated 
for the model. 

1. Skin Friction Drag 

The generic high velocity skin friction drag coefficient used in Chapter II is not 
sufficiently accurate for the performance calculations because of the larger range of 
velocities. The velocity varies from three thousand meter per second to as low as five 
hundred meters per second. The projectile's velocity is still the dominant factor in the 
numerator of the Reynold's number equation, equation (2.24) for the nose body and 
equation (4.9) for the fins. 

Ren/b = pvL T /p (2.24) 
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Re f =pvC ma /jj, (4.9) 

The skin friction drag coefficients calculated with equation (2.25), using the appropriate 
Reynolds number for the nose/body or fins, does not account for compressible flow at 
supersonic speeds. 

C f = .455 [Logi 0 (Re)]‘ 2 ' 58 (2.25) 

The compressible flow factor is calculated with equation (4.10). 

fcf={l + [(Y-l)/2]M 2 }- 467 (4.10) 

Even with the compressible flow factor correction the skin friction coefficient is still 
close to .002. Even a small difference, however, over a long distance, can make a 
difference of kilometers in range, reducing the accuracy of the model. 

The skin friction drag coefficients must be prepared for combination with the 
other drag coefficients for use in the equations of motion by multiplying them by the 
appropriate wetted area and dividing by the reference area, the projectile's maximum 
cross sectional area. These hybrid skin drag coefficients are the normal coefficients 
referenced to a common reference area in order simplify the computer programming and 
to facilitate the notation in this thesis. They are defined in equations (4.11) and (4.12). 

Cfb = Cfl,f c f(S no / s +S b s)/S css (4.11) 

Cff^CfffcfSwf/Scss (4.12) 

2. Form Drag 

The form drag coefficients are determined by the methods from Reference (1). 
The form drag coefficient for the nose/body is calculated with the sane equations as in 
Chapter II. 

C Fb = (Ap/q) {1 - [(196(L n /d) 2 -16)/( 14(M+18)(L n /d) 2 )]} (2.28) 
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Ap/q = [.083+(.096/M 2 )](ct/1 O) 1 ' 69 

(2.29) 

o = tan *(l/2f) 

(2.30) 

The form drag coefficient for the fins is equation (4.13). 


C Ff =(5.33x 2 Sf)/p(S f -.5Seff) 

(4.13) 


The biconvex configuration of the fins accounts for the 5.33. If a different wing 
configuration is desired, e.g. modified double wedge, the correct factor must be 
calculated by the user as shown in Reference (1) and changed in the equations of motion 
as this factor is not in the input section of the program. The terms Sf and S e ff are 
determined using equations (3.3) and (3.21) from Chapter HL 

Sf=(Q-Ct)b + 2Ctb (3.3) 

Seff = S f - {[(1.1452tan(7t/2-A-p))/(tan(7i/2-A-|A)+tan(A-Ji/2-p))] 2 + (3.21) 

[(1.1452tan(A-7i/2-p))/(tan(7c/2-A-p)+tan(A-7c/2-p))] 2 }tan(rc/2-A-p) 

3. Boattail and Base Drag 

The boattail drag, C bt , and base drag, C b , coefficients are empirical parameters 
presented graphically in Reference (1). Approximations are calculated in Appendix (D) 
with the other approximation formulas. The quantities coefficients determined by the 
approximations have to be adjusted to the reference area by dividing by the maximum 
projectile cross sectional area. 

4. Drag Factor 

The drag factor used in the equations of motion is a sum of the drag coefficients 
that have been referenced to the projectile's maximum cross sectional area multiplied by 
the dynamic head and reference area. Equation (4.14) presents the drag factor with the 
dynamic head shown in its component parts. 

CD = (Cft'+Cff+C Fb +C F f+C bt +C b )S C ssPv 2 /2 (4.14) 
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A review of the performance program will reveal that this equation never appears. 
The factors from the drag calculations and equations of motion are at times intermingled 
in the performance program in order to reduce the number of calculations and thus reduce 
the run time of the program. 

D. EQUATIONS OF MOTION 

The equations of motion in the performance program were originally written for 
an assignment in SE 4021 involving the super gun. SE 4021 is a computer simulation 
course in the Physics Department at the Naval Postgraduate School taught at the time by 
Dr. William Colson. The class assignment used a simple model of the atmosphere, a set 
projectile geometry, and a set drag coefficient based on a "bullet" shape. While this 
thesis went into far greater detail actually calculating those quantities, the actual 
equations of motion are not changed from the original assignment. 

1. Velocity 

The initial velocity is inputted by the user as the muzzle velocity. This velocity is 
decomposed into its lateral and vertical components in equations (4.15) and (4.16) by 
using the firing angle specified by the user as the initial angle of the trajectory, <J>. 

v x = vcos<|> (4.15) 

v y = vsin<|) (4.16) 

The performance program calculates the velocity at each time as a change from 
the previous time by a time step. The time step is calculated by dividing the muzzl e 
velocity by the acceleration of gravity and multiplying by a user selected factor as in 
equation (4.17). 

dt = .0001v/g (4-17) 
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Determining the time step in this fashion allows the time step to be of reasonable length 
for the selected muzzle velocity. Using a small factor increases the number of time steps, 
which improves the accuracy of the model, but slows down the run time of the program. 

The lateral and vertical velocities are calculated seperately with equations (4.18) 
and (4.19). 

v x (t) = v x (t-dt) - kv(t-dt)v x (t-dt)dt (4.18) 

v y (t) = Vy(t-dt) - gdt - kv(t-dt)v y (t-dt)dt (4.19) 

The vertical velocity is reduced not only by the effects of drag, but also by gravity. The 
non-subscripted velocity is the total velocity calculated with equation (4.20). 

v = (v x 2 + v y 2 ) 1/2 (4.20) 

The factor k is related to the factor CD, calculated in the previous section, by equation 
(4.21). 

k = CD/mv 2 = (C fb '+C ff '+C F b+CFf+Cbt+Cb)S css p/2m (4.21) 

Using k instead of CD in the equations of motion greatly simplifies the formulas. 
Dividing CD by m converts the drag force into an acceleration. Multiplying an 
acceleration by a time step provides a velocity. Dividing CD by the velocity squared is a 
contrivance to avoid having to track the angle of the trajectory at each time step. In 
eq uatio ns (4.18) and (4.19) the v 2 , which came from the dynamic head calculation, is 
replaced by v(t)v x (t-dt) and v(t)v y (t-dt), as appropriate, which have the angle information 
imb edded in the subscripted velocity as developed for the lateral case in equation (4.22). 

v x = vcos<|> .'. v^osfj) = v(vcos<j>) = w x (4.22) 
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Altitude (km) 


2. Position 


The lateral and vertical positions are once again calculated separately with 
equations (4.23) and (4.24) by adding the change in position since the previous time step 
to the previous time step's position. 

x(t) = x(t-dt) + v x dt (4.23) 

y(t) = y(t-dt) + Vydt (4.24) 

E. RESULTS 



The projectile has a maximum range of one hundred and eighty two kilometers, as 
seen in Figure (16) above, when fired with a fifty seven degree launch angle. That range 
is approximately eight times the maximum range, twenty two kilometers, of a 
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conventional five inch gun. Figure (17) below shows that the projectile only takes four 
minutes and eight seconds to reach a target at its maximum range and eight seconds to 
reach a target at the maximum range of a conventional five inch gun. 


Figure (17) Impact Time 



The kinetic energy of the projectile at impact is also substantially larger than a 
conventional five inch projectile. Figure (18) on the next page shows the energy of the 
projectile at impact for ranges from zero to maximum. There are two energies for almost 
all the ranges because the launch velocity is held constant and the launch angles are 
varied between zero and eighty nine and a half degrees. An interesting feature of Figure 
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(18) is that for ranges greater than sixteen kilometers the projectile has greater kinetic 
energy if shot at the high angle rather than the direct angle. 



Two mechanisms are at work to improve the performance of the high angle shots. 
The first mechanism is that the high angle shots reach the higher atmosphere, which has 
lower density and thus less drag, quicker so they do not slow down as rapidly. The 
second mechanism is that the higher angle of fire shots also descend at a higher angle so 
gravity helps recoop some of the drag loss, i.e. the vertical component of velocity is a 
larger part of the total velocity. 
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Some of the points of interest on Figure (18) are the maximum energy, the energy 
crossover point between direct fire and high angle fire, the minimum energy, and the 
energy at maximum range. At the muzzle of the rail gun the projectile has a maximum 
kinetic energy of three hundred and seven megajoules, this is equivalent to the energy 
released by detonating two hundred and forty four pounds of TNT. The crossover point 
where the kinetic energy of a projectile fired at a shallow angle (direct fire) and a high 
angle occurs at about sisteen kilometers. The kinetic energy at sixteen kilometers is 
twenty seven megajoules, this is equivalent to twenty one pounds of TNT. Referring to 
Figure (17), the projectile takes either eleven or three hundred and twenty five seconds to 
arrive. The minimum kinetic energy of six point eight megajoules, equivalent to only 
five and a half pounds of TNT, occurs at fifty kilometers with a launch angle of twenty 
one and a half degrees. The projectile reaches fifty kilometers in one minute. As 
mentioned above the projectile reaches its maximum range of one hundred and eighty 
two kilo meters in four minutes and eight seconds. The projectile has eighteen 
magajoules of kinetic energy when it impacts at its maximum range, this energy is 
equivalent to fifteen pounds, of TNT. 

The projectile often has the kinetic energy equivalent to an amount of TNT 
substantially larger than the explosive payload of the projectile. The destructive power, 
however, may be substantially less than the TNT because there is no mechanism to 
efficiently transfer the projectile's kinetic energy to the target. 
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V. CONCLUSIONS AND RECOMMENDATIONS 

An aerodynamically stable projectile can be made for a rectangular barrel. Figure 
(1) shows the basic configuration and dimensions of the projectile developed and then 
used as an example in this thesis. The projectile performs very well, but the performance 
model is for the best expected muzzle velocity and with a mass that may be too large for 
practical shipboard deployment. The real benefit of this thesis is the three computer 
programs that will allow future users to quickly analyze different projectile 
configurations to meet evolving requirements. 

A great deal can be done to improve the capabilities of the programs developed in 
this thesis. The development and incorporation of the equations for the subsonic regime 
into the models, with a transsonic approximation to link the subsonic and supersonic 
regimes, will allow a greater variety of materials to be considered. Different fin types, 
i.e. double wedge or modified double wedge, can also be incorporated into the models to 
give users a choice. 

A proper heat transfer study of the projectile should be done. The information 
from a heat transfer study can be used to determine both the fineness ratio at which the 
ogive nose would melt and what materials can be used in the construction of the 
projectile. The form drag of the nose is dependent on the fineness ratio. The 
performance of the projectile can be improved by making the nose pointier to reduce the 
drag, but if the nose is too pointy it may melt. Different materials may be used for 
different portions of the projectile if the temperature of a desired component is known not 
to exceed the melting point of the components constituent material. 
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Wind tunnel testing could verify the drag equations used or provide information 
to improve their accuracy. Hypersonic wind tunnels are very expensive to operate, but 
the proposed projectile flies most of its trajectory at velocities easily achieved in standard 
wind tunnels. A study of the true effects of the flat contact surfaces would determine if 
the assumptions made in calculating the drag of the projectile's body sufficiently 
accounted for the airflow separation and turbulence that the flat contact surfaces will 
probably cause. 

The most difficult study required may well be determining what explosive 
compound to use in the projectile and how it can be fused. The explosive compound 
used will have to be able to undergo the shock of the rapid acceleration experienced in 
the barrel of the rail gun and the high temperatures the projectile may reach without 
detonating. Electrical fuses may be damaged, or worse, detonated by the very high 
currents encountered in the barrel. An impact fuse in the base may work, but a thorough 
study should be done to insure the projectile does not detonate in the barrel of the rail 
gun. 

Experiments with rail guns are on going. This thesis provides tools to help in the 
design of the total weapon system. 
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APPENDIX (A) MATLAB MODEL FOR NOSE/BODY 
CONFIGURATION COMPARISON 


% Inputs 


stin=.5; % skin thickness in inches 

rhos=19200.0; % density of metal in kg/m A 3 (tungsten) 

rhoe=1520. 0; % density of explosive in kg/m A 3 (C4) 

v0=3000; % initial shell velocity in m/s 

cwin=2; % width of contact area in inches 

ew=5.5; % explosive weight in kg 

f=l; % fineness ratio, (nose length)/diameter 

% Constants 

rho=l .03; % ai'r density at sea level in kg/m A 3 

a=346; % speed of sound in air at seal level/standard day in m/s 

% Calculated Initial Values 

M0=v0/a; 

q=.5*rho*(vO A 2); % kg/m-s A 2 

stm=stin*.025416; % skin thickness in meters 

cwm=cwin*.025416; % width of contact area in meters 

% Drag Factors 

Cf=.002; % skin friction drag coefficient from AA3710 notes pg 6-4 
sigmasph=45; % semivertex angle in degree for sphere 

sigmaogi=57.2958*atan(l/(2*f)); % semivertex angle in degrees for ogive 
presssph=(.083+(.096/(M0 A 2)))*((sigmasph/10) A 1.69); 

% from AA3710 notes pg 12-2 
pressogi=(.083+(.096/(M0 A 2)))*((sigmaogi/10) A 1.69); 

CDFsph=presssph*(1-(((392*(.5 A 2))-32)/(26*(MO+18)*(.5 A 2)))); 
CDFogi=pressogi*(1-(((392*(f A 2))-32)/(26*(MO+18)*(f A 2)))); 

% Spherical with Wedge Insert 

srinw=Linspace(2,6,30); % radius in inches 
sdinw=(2.*srinw)+cwin; % Projectile width in inches 
srmw=srinw*.0254166667; % radius in meters 
snvw=(((2/3)*pi).*(srmw. A 3))+(((pi/2)*cwm).*(srmw. A 2)); 

% Nose volume in m A 3 

snww=rhos.*snvw; % Weight of nose in kilograms 
spcsw=(pi.*(srmw. A 2))+((2*cwm).*srmw); 

% Projectile cross section in meters A 2 
spcsiw=spcsw./.000646; % Projectile cross section in in 2 
secsw=(pi.*((srmw-stm). A 2))+((2*cwm).*(srmw-stm)); 

% Explosives cross section in m A 2 
selw=ew./(rhoe.*secsw); % Length of explosive in meters 
spwkw=(rhos.*(((spcsw-secsw) .*selw) + (stm.*spcsw)))+ew+snww; 

% Projectile weight in kilograms 
spwlw=spwkw.*2.2; % Projectile weight in lbm 

stlmw=srmw+selw+stm; % Total length of projectile in meters 
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stliw=stlmw/.0254166667; % Total length in inches 

ssamw"(2.*pi.*(srmw. A 2)) + (pi.*srmw.*cwm)+(2.*pi.*srmw.*(selw+stm) ) ... 

+(2.*(selw+stm).*cwm); % Projectile Surface Area in m A 2 

ssaiw=ssamw./.000646; % Projectile Surface Area in in A 2 
ssfdmw=Cf*ssamw*q; % skin friction drag 
sfdmw==CDFsph*spcsw*q; % form drag 

% Spherical with Shaved Sides 

srins=Linspace(2.0,6,30); % radius in inches 

sdins=(2.*srins); % Projectile width in inches 

srms=srins*.0254166667; % radius in meters 

snvs=(((2/3)*pi).*.(srms. A 3))-((cwm. *(srms. A 2)*pi.*.5)-.*. 

(cwra.*(srms. A 2) .*acos(cwm. / (4.*srms))) - (4.*(srms. A 3).*(sqrt(1-.. . 

(((cwm./(4.*srms)). A 2)))))+(4.*(srms. A 3))-((srms.*cwm./24) .. 

((cwm. /(2.*sin(acos(cwm./(16.*srms))))) + (cwm./(2.*sin(acos(cwm./... 
(8.*srms)))))+(3.*cwm./(2.*sin(acos{(3.*cwm). / (16.*srms)))))+... 
(cwm./(2.*sin(acos(cwm. /(4.*srms))))))) + ((cwm./(96.*srms)) .*... 

(((cwm. A 3)./(128.*sin(acos(cwm./(16.*srms))))) + ((cwm. A 3)./(32.*.. . 

sin(acos(cwm./(8.*srms)))))+((27.*(cwm. A 3))./(128.*sin(acos((3.*cwm)... 
./(16.*srms))))) + ((cwm. A 3)./(8.*sin(acos(cwm./(4.*srms) ))))))); 

% Nose volume in m A 3 

snws=rhos.*snvs; % Weight of nose in kilograms 
spcss—(pi.*(srms. A 2))-(2.*(((srms. A 2).*((pi/2)-... 

acos (.5*cwm./srms) ) ) - ( ( (srms- (srms.Mi- 
sin (acos(.5.*cwm./srms)))))... 

- A 2).*cot(acos(.5*cwm. /srms))))); 

% Projectile cross section in meters A 2 
spcsis=spcss./.000646; % Projectile cross section in in A 2 
secss=(pi.*((srms-stm). A 2))-(2.*((((srms-stm). A 2).*((pi/2)-... 
acos (. 5*cwm. / (srms-stm) )))-((( (srms-stm) - ( (srms-stm) . * . . . 

(1-sin(acos(.5.*cwm./(srms-stm) ))))). A 2).*. 

cot(acos(.5*cwm./(srms-stm) ))))); % Explosives cross section in m A 2 

sels=ew./(rhoe.*secss); % Length of explosive in meters 
spwks=(rhos.*(((spcss-secss).*sels)+(stm.*spcss)))+ew+snws; 

% Projectile weight in kilograms 
spwls=spwks.*2.2; % Projectile weight in lbm 
stlms=srms+sels+stm; % Total length of projectile in meters 
stlis=stlms/.0254166667; % Total length in inches 

ssams=(2.*pi.*(srms. A 2))-(2.*pi.*(((-(srms. A 2). 12 ) .*log((.5.*cwm) +_ 

sqrt((srms. A 2)+ (.25.*(cwm A 2))))) + (.25.*cwm.*sqrt((srms. A 2) + (.25.*. . . 
(cwm A 2))))+<.5.*(srms. A 2).*log(srms))))+(pi.*((cwm.*.5). A 2)) + ... 

(4.*srms.*(sels+stm).*asin(.5.*cwm./srms)) + (2.*cwm.*(sels+stm) ) ; 

% Projectile Surface Area in m A 2 
ssais=ssams./.000646; % Projectile Surface Area in in A 2 
ssfdms=Cf*ssams*q; % skin friction drag 
sfdms=CDFsph*spcss*q; % form drag 


% Ogival with Wedge Insert 

orinw=Linspace(2,6,30); % radius in inches 
odiriw=(2.*orinw)+cwin; % Projectile width in inches 
ormw=orinw.*.0254166667; % radius in meters 
onlw=2.*ormw.*f; % nose length in meters 
oRmw=((onlw. A 2)+(ormw. A 2))./(2.*ormw); 
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% Radius of curvature of nose in meters 
thetaw=acos(onlw./oRmw); 

onvw=(pi-*((onlw.*((2.*(oRmw. A 2))-(2.*oRmw.*ormw)+(ormw. A 2)))-... 

((onlw. A 3)./3)-(onlw.*(oRmw-ormw).*sqrt((oRmw. A 2)-(onlw. A 2)>)••• 

-{(oRmw-ormw).*(oRmw. A 2).*asin(onlw./oRmw))))+(cwm.*(((oRmw. A 2)... 
.*((pi/2)-thetaw))-(((oRmw-ormw). A 2).*cot(thetaw)))); 

% Nose volume in m A 3 

onww=rhos.*onvw; % Weight of nose in kilograms 
opcsw=(pi.*(ormw. A 2))+(2.*cwm.*ormw); 

% Projectile cross section in meters A 2 
opcsiw=opcsw./.000646; % Projectile cross section in in 2 
oecsw=(pi.*((ormw-stm). A 2))+(2.*cwm.*(ormw-stm)); 

% Explosives cross section in m A 2 
oelw=ew./(rhoe.*oecsw); % Length of explosive in meters 
opwkw=(rhos.*(((opcsw-oecsw).*oelw)+(stm.*opcsw)))+ew+onww; 

% Projectile weight in kilograms 
opwlw=opwkw.*2.2; % Projectile weight in lbm 

otlmw=(f.*2.*ormw)+oelw+stm; % Total length of projectile in meters 

otliw=otlmw/.0254166667; % Total length in inches 

osamw=(2.*pi.*((2.*f.*ormw.*oRmw)-(oRmw.*(oRmw-ormw).*asin... 

(2.*f.*ormw./oRmw))))+(2.*cwm.*oRmw.*((pi/2)-thetaw))+... 

(2.*pi.*ormw.*(oelw+stm)) + (2.*cwm.*(oelw+stm)); 

% Projectile Surface Area in m A 2 
osaiw=osamw./.000646; % Projectile Surface Area in in A 2 
osfdmw=Cf*osamw*q; % skin friction drag 
ofdmw=CDFogi*opcsw*q; % form drag 


% Ogival with Shaved Sides 

orins=Linspace(2.0,6/30); % radius in inches 
odins=(2.*orins); % Projectile width in inches 
orms=orins*.0254166667; % radius in meters 
onls=2.*orms.*f; % nose length in meters 
oRms=((onls. A 2)+(orms. A 2))./(2.*orms); 

% Radius of curvature of nose in meters 
thetas=acos(onls./oRms); 

onvs=(pi.*((onls.*((2.* (oRms. A 2)) - (2. *oRms. *orms) + (orms. A 2)))-... 

((onls. A 3)./3)-(onls.*(oRms-orms).*sqrt((oRms. A 2)-(onls. A 2)))... 

-((oRms-orms).*(oRms. A 2).*asin(onls./oRms))))-.((cwm.*(oRms. A 2)*... 
pi.*.5)-(cwm.*(oRms. A 2).*acos((f.*cwm)./(2.*oRms)))-((2./f).*... 
(oRms. A 3) .* (sqrt (1- (((f.*cwm) ./ (2. *oRms)) . A 2))) ) + ((2.* (oRms. A .. . 

3)) ./f)-((f.*oRms.*cwm./12).*((cwm./(2.*sin(acos((f.*cwm) ./(8.*... 
oRms)))))+(cwm./(2.*sin(acos((f.*cwm)./(4.*oRms)))))+(3.*cwm./... 

(2.*sin(acos((3.*f.*cwm)./(8.*oRms))))) + (cwm. /(2.*sin(acos((f.*... 
cwm) ./(2.*oRms))))))) + (((cwm.*(f. A 3))./(12.*oRms)).*(((cwm. A 3)./... 
(128.*sin(acos((f.*cwm)./(8.*oRms)))))+((cwm. A 3)./(32.*sin(acos... 

((f.*cwm) ./(4.*oRms))))) + ((27.*(cwm. A 3))./(128.*sin(acos{(3.*f.*..• 

cwm)./(8.*oRms))))) + ((cwm. A 3)./(8.*sin(acos((f.*cwm) ./(2.*oRms)))))))) 
% Nose volume in m A 3 

onws=rhos.*onvs; % Weight of nose in kilograms 

opcss=(pi.*(orms. A 2))-(2.*(((orms. A 2).*((pi/2)-acos(.5*cwm. /orms)))-.. 
(((orms-(orms.*(1-sin(acos(.5.*cwm./orms))))). A 2).*cot(acos(.5.*... 
cwm./orms))))); 

% Projectile cross section in meters A 2 
opcsis=opcss./.000646; % Projectile cross section in in A 2 
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oecss=(pi.*((orms-stm). A 2 ))- ( 2 .* ( ( ( (orms-stm) . A 2 ).*((pi/ 2 ) 
acos (.5*cwm. / (orms-stm) )))-((((orms-stm) - ((orms-stm) . * . .. 

( 1 -sin(acos(. 5 .*cwm./ (orms-stm)))))). A 2 ) .*... 

cot(acos(.5*cwm./(orms-stm)))))); % Explosives cross section in m A 2 

oels=ew./(rhoe.*oecss); % Length of explosive in meters 
opwks=(rhos.*(((opcss-oecss) . *oels)+ (stm. *opcss)))tew+onws; 

% Projectile weight in kilograms 
opwls=opwks.*2.2; % Projectile weight in lbm 

otlms-(f .*2.*orms) -f-oels+stm; % Total length of projectile in meters 
otlis=otlms/.0254166667; % Total length in inches 
osams=(2.*pi.* ( (2. *f. *orms. *oRms) 7 (oRms. * (oRms-orms) ,.*asin. . . 

(2. *f. *orms. /oRms) ) ) ) - (oRms. * ( (pi. *cwm) - (2 . *cwm. *acos ( (f. *cwm) . / .. . 

(2.*oRms)))-((4.*oRms.*sqrt(l-(((f.*cwm) ./(2.*oRms)). A 2)))./f) +_ 

(4.*oRms./f)))+(2.*({(oRms. A 2).*((pi./2)-acos(cwm./oRms) ))-... 

( ( (oRms- (cwm/2) ) . A 2) . *cot (acos (cwm. /oRms))))) + (4. *orms. * (oels+. .. 
stm) . *asin (. 5. *cwm. /orms) ) + (2. *cwm. * (oels+stm)) ; 

% Projectile Surface Area in m A 2 
osais=osams./.000646; % Projectile Surface Area in in A 2 
osfdms=Cf*osams*q; % skin friction drag 
ofdms=CDFogi*opcss*q; % form drag 

figure(1) 

plot(sdinw,spcsiw, T k- 1 ,sdins,spcsis, ! k. 1 ),grid 
title(’Cross Section vs Width*) 

xlabel( T Projectile Width (in) r ),ylabel('Projectile Cross Section 
(in A 2)*) 

text(6.1,110,* Solid = w/ Insert 1 ) 
text(6.1,105, T Dotted ~ w/ Shaved Sides’) 
figure(2) 

plot(sdinw,spwlw, 1 k- f ,sdins,spwls,’k. 1 ,odinw,opwlw,* k.-. f ,odins,... 

opwls, r k—’),grid 
title{’Weight vs Width’) 

xlabel( 1 Projectile Width (in)’), ylabel(’Projectile Weight (lbm)’) 

text (6.1, 654, ’ Solid =* Spherical w/ Insert’) 

text(6.1,614, ’Dotted = Spherical w/ Shaved Sides’) 

text(6.1,574,’Dash/Dot = Ogival w/ Insert’) 

text(6.1,534,’Dashed = Ogival w/ Shaved Sides’) 

figure(3) 

plot (sdinw, stliw, ' k- *, sdins, stlis, ’ k. 1 , odinw, otliw, ’ k. -. ’, odins, . .. 

otlis, 1 k— 1 ),grid 
title(’Length vs Width’) 

xlabel(’Projectile Width (in) 1 ), ylabel(’Total Length (in)’) 

text(6.1,41,’Solid = Spherical w/ Insert') 

text(6.1,38.5,’Dotted = Spherical w/ Shaved Sides’) 

text(6.1,36,’Dash/Dot = Ogival w/ Insert’) 

text(6.1, 33.5, ’Dashed = Ogival w/ Shaved Sides’) 

figure(4) 

plot(spwlw,stliw,’k-’,spwls,stlis,’k.’,opwlw,otliw,’k.,opwls,... 

otlis,’k—’),grid 
title(’Length vs Weight’) 

xlabel(’Projectile Weight (lbm)’), ylabel(’Total Length (in)’) 

text(205,41,'Solid = Spherical w/ Insert’) 

text(205,38.5, ’Dotted = Spherical w/ Shaved Sides’) 

text(205,36,’Dash/Dot = Ogival w/ Insert') 

text(205, 33.5,’Dashed = Ogival w/ Shaved Sides’) 

figure(5) 

plot (spwlw, spcsiw, ’ k- ’, spwls, spcsis, ’ k. ’, opwlw, opcsiw, ’ k.-. ’, opwls, . . . 
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opcsis,'k—'),grid 
title('Cross Section vs Weight') 

xlabelCProjectile Weight (Ibm)'),ylabel('Projectile Cross Section 
(in^2)’) 

text(215,37,'Solid = Spherical w/ Insert’) 
text(215,30,'Dotted = Spherical w/ Shaved Sides') 
text(215,23,'Dash/Dot = Ogival w/ Insert') 
text(215,16,'Dashed = Ogival w/ Shaved Sides') 
figure(6) 

plot(spcsiw,ssaiw,'k-',spcsis,ssais,'k.',opcsiw,osaiw,'k.',opcsis,... 
osais,'k—'),grid 

title('Surface Area vs Cross Section') 
xlabelCProjectile Cross Section (in-2)') 
ylabel('Projectile Surface Area (in-2)') 
text(21, 480,'Solid = Spherical w/ Insert') 
text(21,460,'Dotted = Spherical w/ Shaved Sides') 
text(21, 440,’Dash/Dot = Ogival w/ Insert') 
text(21,420,'Dashed = Ogival w/ Shaved Sides') 
figure(7) 

plot(ssfdmw,sfdmw./le5,'k- 

’,ssfdms,sfdms./Ie5,'k.',osfdmw,ofdmw./Ie5,... 

' k.-. ', osfdms, ofdms. /Ie5, ' k—'), grid 
title('Form Drag vs Skin Friction Drag at Sea Level, V0=3000 km/s') 
xlabel('Skin Friction Drag (kg-m/s A 2)'),ylabel('Form Drag x 100,000 
(kg-m/s-2)') 

text(1750,1.56,'Solid = Spherical w/ Insert') 
text(1750,1.41,'Dotted = Spherical w/ Shaved Sides') 
text(1750,1.26,'Dash/Dot = Ogival w/ Insert') 
text(1750,1.11,'Dashed = Ogival w/ Shaved Sides') 
figure(8) 

plot(sdinw,ssfdmw,'k-',sdins,ssfdms,'k.',odinw,osfdmw,'k.-.',odins,... 
os f dms,'k—'),grid 

title('Skin Friction Drag vs Width at Sea Level, V0=3000 km/s') 

xlabel('Projectile Width (in)'),ylabel('Skin Friction Drag (kg-m/s A 2)') 

text(6.1,2900,'Solid = Spherical w/ Insert') 

text(6.1,2800,'Dotted = Spherical w/ Shaved Sides') 

text(6.1,2700,'Dash/Dot = Ogival w/ Insert') 

text(6.1,2600,'Dashed = Ogival w/ Shaved Sides') 

figure(9) 

plot(sdinw,sfdmw./Ie5,'k-',sdins,sfdms./Ie5,'k.',odinw,ofdmw./Ie5,... 

'k.-.',odins,ofdms./Ie5,'k—'),grid 
title(Form Drag vs Width at Sea Level, V0=3000 km/s') 

xlabel('Projectile Width (in)'),ylabel('Form Drag/100,000 (kg-m/s A 2)') 

text(5,2.41,'Solid = Spherical w/ Insert') 

text(5,2.24,'Dotted = Spherical w/ Shaved Sides') 

text(5,2.07,'Dash/Dot = Ogival w/ Insert') 

text(5,1.90,'Dashed = Ogival w/ Shaved Sides') 


If a color printer is available different colored lines are a better presentation on the plots. 
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APPENDIX (B) PERFORMANCE PROGRAM (in C) 

/*Performance*/ 

#include <stdio.h> 

#include <math.h> 

main(){ , 

float g,C, rhoo, rho. A, v, vx, vy,x,y, theta, t,dt,k, h,xf, yf,pi,r,f; 
float a,M,T,Reb,Ref,CDFb,Cfb,CDfb,pr,gamma,R,mu,SA,pb,dc,dee,er; 
float CDbt,CDb,nb,lb,Sb,rin,stin,rhom,rhoe,stm,nl,nv,nw,br,lbt; 
float ecs,el,ew,nlro,nlm2,nlm3,tlm, 1, sigma, ein,r2,r3, sigmac, ewin; 
float cr,ct,lamdam,mac,tau,tmac,rf,phi,bin,b,fv,fw,fsa,m, cw3; 
float ibrh,ibrw,omega,pf,Cff,CDff,CDFf, beta,Mu, seff,sf,et,nbt; 
float lamda,psi, Psi, Dl, D2, D3, comf, rn, m2, m3, thetan, crp, f3, cw; 
int j; 

/* INPUTS */ 

rin=2.5; /* radius in inches */ 
stin=.5; /* skin thickness in inches */ 

' cwin=2; /* width of contact area in inches */ 
rhom=19200.0; 

/* density of metal in kg/m A 3 (tungsten in example) */ 
rhoe=1520.0; 

/* density of explosive in kg/m A 3 (04 in this example) */ 
ew=5.5; /* desired explosive payload in kg */ 
v=3000.0; /* muzzle velocity in m/s */ 

theta=.986111; /* launch angle in radians (56.5 degrees) */ 
tau^.03'; /* fin thickness to chord ratio */ 
crp=-6; /* fin percentage of total body length */ 
lamda=.585686; /* leading edge sweep angle */ 

f=l; /* nose fineness ratio, (nose length)/(body diameter) */ 

/* CONSTANTS */ 

rhoo=1.03; /* air density at sea level in kg/m A 3 */ 
pi=3.141593; 

g=9.81; /* gravity in m/s A 2 */ 
gamma=1.4; /* gamma for air */ 

R=287; /* gas constant for air in SI units*/ 

- mu=1.8e-5; /* coefficient of viscosity */ 
a=340.17; 

/* Speed of sound at sea level on standard day in m/s */ 

/* INITIAL VALUES */ 
x=0.0; /* initial position */ 
y=0.0; /* initial altitude */ 
t=0.0; /* initial time */ 
j=0; /* counter variable */ 

/* CALCULATED INITITIAL VALUES AND CONSTANTS */ 
vx=v*cos(theta); /* initial lateral velocity in m/s */ 
vy=v*sin(theta); /* initial vertical velocity in m/s */ 

M=v/a; /* Mach Number for muzzle velocity at sea level */ 

/* Projectile Dimensions and Parameters */ 
r=rin*.025417; /* radius in meters */ 
r2=r*r; 
r3=r*r*r; 
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stm=stin*.025417; /* skin thickness in meters */ 
cw=cwin*.025417; /* width of contact area in meters */ 
cw3=cw*cw*cw; 

nl=f* 2 *rin; /* nose length in inches */ 
nlm=f* 2 *r; /* nose length in m */ 
nlm 2 =nlm*nlm; 
nlm3=nlm*nlm*nlm; 

rn=((nlm*nlm)+r 2 )/( 2 *r); /*radius of curviture of nose in m */ 

rn 2 =rn*rn; 

rn3=rn*rn*rn; 

thetan=acos(nlm/rn); /* angle of arc swept by nose */ 
ein=rin-stin; /* explosive radius in inches */ 
er=ein*.025417; '/* explosive radius in meters */ 
f3=f*f*f; 

nv=(pi*((nlm*((2*rn2)-(2*rn*r)+r2))-(nlm3/3)-((rn-r)*nlm* 

sqrt(rn 2 -nlm 2 ))-((rn-r)*rn 2 *asin(nlm/rn))))-((cw*rn 2 * 
pi*. 5 )-(ew*rn 2 *aeos((f*cw)/( 2 *rn)))-(( 2 /f)*rn 3 *sqrt( 1 - 

(((f*cw)/(2*rn))*((f*cw)/(2*rn)))))+((2*rn3) /£)-{ (f*rn* 
cw/ 12 )*((cw/( 2 *sin(acos((f*cw)/( 8 *rn))))) + (cw/( 2 *sin 
(acos((f*cw)/(4*rn)))))+(3*cw/( 2 *sin(acos((3*f*cw)/ 

( 8 *rn)))))+(cw/( 2 *sin(acos((f*cw)/( 2 *rn)))))))+(((cw*f3)/ 

(12*rn))*((cw3/(128*sin(acos((f*cw)/(8*rn)))))+(cw3/ 

(32*sin(acos((f*cw)/(4*rn)))))+((27*cw3)/(128*sin(acos 
((3*f*cw)/( 8 *rn)))))+(cw3/( 8 *sin(acos((f*cw)/(2*rn)))))))); 

/* nose volume in m A 3 */ 
nw=rhom*nv; /* nose mass in kilograms */ 

A=(pi*r2)-(2*((r2*((pi/2).-acos(.5*cw/r)))-(((r-(r*(1-sin(acos 
. (.5*cw/r)))))*(r-(r*( 1 -sin(acos(.5*cw/r))))))*(.5*cw/r)* 
( 1 /sin(acos(.5*cw/r)))))); 

/* projectile cross section in m A 2 */ 
ecs=(pi*((r-stm)*(r-stm)))-( 2 *(( ((r-stm) *(r-stm))*((pi/ 2 )-acos 
(. 5 *cw/(r-stm))))-((((r-stm) -((r-stm) *( 1 -sin(acos(. 5 *cw/ 
(r-stm))))))*((r-stm) -((r-stm) *( 1 -sin(acos(. 5 *cw/ 

(rstm)))))))*(.5*cw/(r-stm))* (1/sin(acos(.5*cw/(r- 
stm))))))); 

/* explosives cross section in m A 2 */ 
el=ew/(rhoe*ecs); /* length of explosive in m */ 
tlm=nlm+el+stm; /* total length of- projectile in m */ 

1= (5*r* (l-sin'(acos (2/(5*rin))) ) ) ; /* intermediate calculation */ 
SA=(2*pi*((2*f*r*rn)-(rn*(rn-r)*asin(2*f*r/rn))))-(rn*((pi*cw)- 
(2*cw*acos((f*cw)/(2*rn)))-((4*rn*sqrt(1-(( (f*cw)/(2*rn)) 

* ( (f*cw) / (2*rn)))) ) /f) + (4*rn/f) ) ) + (2* ( (m2* ( (pi/2)-acos 
(cw/rn) )) - (( (rn- (cw/2)) * (rn- (cw/2))) * (cw/rn) * (1/sin (acos 
(cw/rn)))))) + (4*r*(el+stm)tasin(.5*cw/r)) + (2*cw*(el+stm)); 
/* projectile surface area in m A 2 */ 
sigma=(360*atan(rin/nl))/pi; 

/* semivertex angle at ogive point in degrees */ 
br=r-(tlm*.014054); 

/* base r in m w/ 8 deg boatail .1 the length of projectile */ 
omega=(pi/ 2 )-asin(cw/( 2 *r)); 

/* half angle of base's curved section */ 

ibrh=br*cos(omega); /* intermediate step */ 
ibrw=br*cos((pi/ 2 )-omega); /* intermediate step */ 

Sb= (2*((br*br*((pi/2)-omega))-(ibrh*ibrh*(cos(omega) 

/sin(omega)))))+( 2 *cw*ibrw); /* base area in m A 2 */ 
cr=crp*tlm; /* root chord in meters */ 

psi=(pi/ 2 )-lamda; /* "inside" sweep angle of leading edge */ 
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b=2*((sqrt(2*((r*sin(acos(.5*cw/r)))*(r*sin(acos 

(.5*cw/r))))))-r); /* span of single fin in meters */ 

ct=(cr-(.5*b/tan(psi))); /* tip chord in meters */ 

Psi=pi-psi; /* "inside sweep 1 angle of tip */ 

sf=(.5*b*(cr-ct))+(ct*b); /* planform area of fin pair in m A 2 */ 

lamdam=ct/cr; /* intermediate calculation */ 

mac= (. 666667*cr* (l+lamdam+ (lamdam*lamdam))) / (1+lamdam); 

/* mean aerodynamic chord in meters */ 
tmac=mac*tau; /* fin thickness at MAC in meters */ 
rf=mac*(l+(tau*tau))/(4*tau); 

/* radius of curviture of biconvex fins at MAC in meters */ 
phi=(pi/2)-asin(mac/(2*rf)); 

/* lower limit of integration for quarter of fin volume */ 
fv=2*b*(((rf*rf)*((pi/2)-phi))-(((rf*rf)-(rf*tmac)+(.25*tmac 
*tmac))*cos(phi)/sin(phi))); 

/* volume of single fin in m A 3 */ 
fw=2*rhom*fv; /* mass of all four fins in kg */ 
fsa=16*rf*((pi/2)-phi)*b; /* surface area of four fins */ 
m=(rhom*(((A-ecs)*el)+(stm*A)))+ew+nw+fw 
/* projectile mass in kg */ 
dt=.0001*(v/g); /* time step for integration */ 

printf("%f, %f, %f,%f,%f\n",x/1000,y/1000,v,t,M); 
while(y>=0 && y<11000){ 

T=288-((6.5*y)/1000);/* Adiabatic Cooling w/ Alt Increase*/ 

a=sqrt(gamma*R*T); 

dc=((y*g)/(R*(T-288.16)))+l; 

dce=pow((T/288.16) f -dc) ; 

rho=rhoo*dce; 

v=sqrt(vx*vx+vy*vy); 

M=v/a; 

beta=sqrt((M*M)-1); 

Mu=asin(l/M); 

Dl=tan(psi-Mu);D2=tan(Psi-Mu);D3=D1+D2; 

/* intermediate calcs */ 

Reb=(rho*v*tlm)/mu; /* Reynolds number for body */ 

Ref=(rho *v*mac)/mu; /* Reynolds number for fins */ 
pb=pow((log(Reb)/log(10)),-2.58); 
pf=pow((log(Ref)/log(10)),-2.58); 

Cfb=.455*pb; /*Yellow book*/ 

Cff=.455*pf; 

comf=pow((1+(((gamma-1)/2)*M*M)),-.467); 

/* compressible flow factor */ 

CDfb=(Cfb*SA*comf)/A; 

CDff=(Cff*fsa*comf)/A; 
sigmac=pow((sigma/10),1.69); 
pr=sigmac*(.083+(.096/(M*M))); 

CDFb=pr*(1-(((196*f*f)-16)/(14*(M+18)*(f*f)))); 
seff=sf-(((1.1452*D1/D3)*(1.1452*D1/D3)*D2)+ 

((1.1452*D2/D3)*(1.1452*D2/D3)*D1)); 

CDFf=(5.33*tau*tau*(sf-(.5*seff)))/(beta*sf) ; 
lbt=500000*((.0555*M*M)+(.006*M)+(-.00255*M*M*M)+.85); 
nbt=((-485*M)*(-485*M) )-(lbt); 

CDbt=(10*((485*M)-sqrt(nbt))/250000)-.00375; 
lb=4000*((.007*M*M) + {.00007*M) + (-.0009*M*M*M)+.85) ; 
nb=((-445*M)*(-445*M))-(lb); 

CDb=110*Sb*((445*M)-sqrt(nb))/(2000*A) ; 
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C=CDfb+CDff+CDFb+CDFf+CDbt+CDb; 
k=(C*rho*A)/(2*m); 
vx=vx-(k*v*vx*dt); 
vy=vy-g*dt-(k*v*vy*dt); 
x=x+vx*dt;y=y+vy*dt;t=t+dt; 

j=j+i; 
if(j>=100){ 

printf("%f,%f,%f,%f,%f\n",x/1000,y/1000,v,t,M); 
j=0;}} 

while(y>=11000 && y<25000){ 

Isothermal Region 
Equations of Motion Loop} 
while(y>=25000 && y<47000){ 

T=216.5+(3*(y-25000)/1000); /* Adiab. Heating w/ Alt Inc* 
Equations of Motion Loop} 
while(y>=47000 && y<53000){ 

Isothermal Region 
Equations of Motion Loop} 
while(y>=53000 && y<79000){ 

T=282.5-(4.5*(y-53000)/1000); /* Adiab. Cooling w/ Alt 

Inc*/ 

Equations of Motion Loop} 
while(y>=79000 && y<=90000){ 

Isothermal Region 
Equations of Motion Loop} 
while(y>=90000 && y<=105000){ 

T=165.5+(4*(y-90000)/10Q0); /* Adiab. Heating w/ Alt Inc*/ 
Equations of Motion Loop} 

while(y>=105000){Drag is Effectively Zero Above this Altitude so 
all the drag calculations are removed from the loop and the 
drag terms are removed from the equations of motion} 
while(y>=90000){ 

Isothermal Region 

T=165.5+(4*(y-90000)/1000); /* Adiab. Cooling w/ Alt Dec*/ 
Equations of Motion Loop} 
while(y>~79000){ 

Isothermal Region 
Equations of Motion Loop} 

' while(y>=53000){ 

T=282.5-(4.5*(y-53000)/1000);/* Adiab. Heating w/ Alt Dec*/ 
Equations of Motion Loop} 
while(y>=47000){ 

Isothermal Region 
Equations of Motion Loop} 
while(y>=25000){ 

T=216.5+(3*<y-25000)/1000); /* Adiab. Cooling w/ Alt Dec*/ 
Equations of Motion Loop} 
while(y>=11000){ 

Isothermal Region 
Equations of Motion Loop} 
while (y>= 5 0) { 

T=288-((6.5*y)/1000); /* Adiabatic Heating w/ Alt Dec*/ 
Equations of Motion Loop} 
xf—X“((vx*y)/vy); 
yf=y-((vy*y)/vy); 

printf (”%f,%f f %f,%f, %f \n f, f xf/1000, yf/1000, v, t f M);} 
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APPENDIX (C) MATLAB STABILITY PROGRAM 

% Stability Study 
% Inputs 

r= 2 . 5 ; % projectile radius in inches 
h= 2 ; % width of contact area in inches 
f=l; % nose fineness ratio 
st=. 5 ; % skin thickness in inches 

rhos=19200.0; % density of metal in kg/m A 3 (tungsten) 

rhoe=1520.0; % density of explosive in kg/m A 3 (C4) 

ew=5.5; % explosive weight in kg 

tau=0.03; % fin thickness ratio 
z=. 6 ; % root chord/total length* 

% Constants 
gamma=l'. 40; 

rho=.00004; % air density in lbm/cubic inch 

T=520; % Standard temperture at sea level in Rankine 

a=(sqrt(gamma*32.2*53.3*T))*12; % Speed of sound in inches/second 

% Nose/Body 

d=(2*r); % Projectile width in inches 

rm=r*.0254166667; % radius in meters 

hm=h*.0254166667; % width of contact areas in meters 

stm=st*.0254166667; % skin thickness in meters 

nl= 2 *rm*f; % nose length in meters 

nli= 2 *r*f; % nose length in inches 

Rm=< (nl A 2)+ (rm A 2)) /(2*rm) ; 

% Radius of curvature of nose in meters 
thetas=acos(nl/Rm); 

nv=(pi*((nl*((2*(Rm A 2))-(2*Rm*rm) + (rm A 2)))-((nl A 3)/3)-(nl*(Rm-rm)*.. 
sqrt((Rm A 2)-(nl A 2)))-((Rm-rm) *(Rm. A 2)*asin(nl/Rm))))-((hm*(Rm A 2)* 
pi*.5)-(hm*(Rm A 2)*acos((f*hm) / (2*Rm)))-((2/f)*(Rm A 3)*(sqrt(1-... 

(((f*hm)/(2*Rm)) A 2))))+((2*(Rm A 3) )/£)-( (f*Rm*hm/12)*((hm/(2*sin.. 
(acos((f*hm)/(8*Rm))))) + (hm/(2*sin(acos((f*hm)/(4*Rm))))) + (3*hm./ 

(2*sin(acos((3*f*hm)/( 8 *Rm))))) + (hm/(2*sin(acos((f*hm)/(2*Rm))))))) + 
(((hm*(f A 3))/(12*Rm))*(((hm A 3)/(128*sin(acos((f*hm)/(8*Rm)))))+.. 
((hm A 3)/(32*sin(acos((f*hm) /(4*Rm))))) + ((27*(hm A 3))/(128*sin(acos 
((3*f*hm)/( 8 *Rm)))))+((hm A 3)/( 8 *sin(acos((f*hm)/(2*Rm)))))))); 

% Nose volume in m A 3 

nw=rhos.*nv; % Weight of nose in kilograms 

pcs=(pi*(rm A 2 ))-( 2 *(((rm A 2 )*((pi/ 2 )-acos(.5*hm/rm)))-(((rm-(rm*( 1 - 
sin... 

(acos(.5*hm/rm))))) A 2)*cot(acos(.5*hm/rm))))) ; 

% Projectile cross section in meters A 2 
pcsi=pcs/.000646; % Projectile cross section in in A 2 
ecs=(pi*((rm-stm). A 2 ))-( 2 *((((rm-stm) A 2 )*((pi/ 2 )-acos(.5*hm/(rm- 
stm))))-... 

((((rm-stm)-((rm-stm)*( 1 -sin(acos(.5*hm/(rm-stm)))))) A 2 )*cot(acos 
(.5*hm/ (rm-stm)))))); % Explosives cross section in m A 2 

% Explosives cross section in m A 2 
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el=ew/(rhoe*ecs); % Length of explosive in meters 
eli=el/.025416667; % Length of explosive in inches 
pwk=(rhos*(((pcs-ecs)*el)+(stm*pcs)))+ew+nw; 

% Projectile weight in kilograms 
pwl=pwk*2.2; % Projectile weight in lbm 

tlm-(f *2*rm)+el+stm; % Total length of projectile in meters 
tli=tlm/.0254166667; % Total length in inches 

sam= (2*pi* ( (2*f *rm*Rm) - (Rm* (Rm-rm) *asin (2*f *rm/Rm)) ) ) - (Rm* ( (pi*hm) - 
(2*... 

hm*acos ( (f *hm) / (2*Rm) ) ) - ( (4*Rm*sqrt (1- ( ( (f*hm) / (2*Rm) ) A 2) ) ) /f) +.. . 
(4*Rm/ f)))+(2*((<Rm A 2)*((pi/2)-acos(hm/Rm)))-(((Rm-(hm/2)) A 2)*cot.., 
(acos(hm/Rm)))))+(4*rm*(el+stm)*asin(.5*hm/rm))+(2*hm*(el+stm)); 

% Projectile Surface Area in m A 2 
sai=sam./.000646; % Projectile Surface Area in in A 2 

%fins 

b=2*((sqrt(2*((r*sin(acos(.5*h/r))) A 2)))-r); % span in inches w/o body 
bm=b*.025416667; % span in meters 
cr=z*tli; % centerline chord in inches 
lamda=.585686; 

% sweep ang.of 34 deg(in rads); keeps LE supersonic down to Mach 1 
theta=(pi/2)-lamda; % "inside" sweep angle of 56 degrees 
ct“(cr-(.5*b/tan(theta))); % tip chord in inches 

Sf=(.5*b*(cr-ct))+(ct*b); % planform area of fin pair in square 

inches 

Sd=(cr A 2)*tan(theta); % area if delta where whole in sqare inches 

Sl=Sf/Sd; * 

S2=Sf/pcsi; 

mac=(2*cr*(l+(ct/cr)+((ct/cr) A 2 )))/(3*(1 +(ct/cr))); 

% mean aerodynamic chord of fins in inches 

macm=mac*.025416667; % mean aerodynamic chord of fins in meters 
Rmf=(((.5*mac) A 2)+((.5*tau*mac) A 2))/(tau*mac); 

% radius of curviture of biconvex fin at mac in inches 
Rmfm=Rmf*.0254166667; 

% radius of curviture of biconvex fin at mac in meters 
thetaf=acos(.5*mac/Rmf); 

fw=2*rhos*bm*(((Rmfm A 2)*((pi/2)-thetaf))-(((Rmfm-(.5*tau*macm)) A 2)... 

*cot(thetaf))); % Weight of single fin in kg 
phi=pi-theta; % "inside” tip angle of 129 degrees 
AR=b A 2/Sf; 

Kl=5.33; 

Apri=2*tau/3; 

% Center of Gravity 

tw=pwk+(2^fw); % weight of entire projectile including fins in kg 
CGn=.707*nli*nw; % CG point for the nose 

CGb=((.5*eli)+nli)*((rhos*(pcs-ecs)*el)+ew); % CG point for body 
CGep=((•5*st)+eli+nli)*(stm*pcs*rhos); % CG point for the end plate 
CGf=((mac*.5)+((l-(z+,l))*tli))*2*fw; 

CG=(CGn+CGb+CGep+CGf)/tw; % CG point of projectile 

% Interference Factors 
rs=r/(r+(b/2)); 

Kbf=(.869*((rs+.65) A 2))-.367; 

Rfb=(.278*((rs+1.3) A 2))+.529; 
kbf=rs; 

kfb= (. 2* ( (rs-. 5) A 2) ) + . 95; 
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M=linspace(1.3,9,250); 
beta=sqrt (M. A 2 -l) ; 
mu=asin(1./M); 
if (beta./f)<=l 

CNan=(-.28.*(((beta./f)-1.9). A 2))+3.58; 

% Approximation over Mach range for chart pg 13-4 
Xn=CG-(tli.*((.125.*beta./f)+.15)); 

% From graph pg 14-8 adjusted to CG 
else 

CNan=(-3.4.*(((f./beta)-.9). A 2))+3.4; 

Xn=CG-(tli.*((-.125.*f./beta)+.425)); 
end 

Dl=tan(theta-mu); 

D2=tan(phi-mu); 

D3=D1+D2^ 

if ((b/ 2 )-((cr+ct).*tan(mu )))<=0 

Seff=Sf-((((1.1452.*D1./D3). A 2).*D2)+(((1.1452.*D2./D3). A 2).*D1)); 
% Effected tip area (split into two terms to account for overlap) 
©Is© 

Seff=((cr A 2).*tan(mu))+((ct A 2).*tan(mu)); 

end 

CNa2df=(4./beta); 

% CNalpha before tip effects 
CNaf=CNa2df-((2.*Seff)./(beta.*Sf)); 

CNa=CNan+(CNaf.*(Kbf+Kfb).*S2) 

CNd=CNaf .* (.(kbf+kfb) *S2) ; 

Xf=-(cr*. 666666 *(SI-(.5.*Seff./Sd)))-(((l-(z+.l))*tli)-CG); 

% CP of fin adjusted for missing area and to CG 

Xbf=-(cr.*((.9888;*beta.*AR)+.35))-(((l-(z+.l))*tli)-CG); 

% From NACA figure 15i and adjusted to CG 

CMa=((CNan.*Xn)+(CNaf.*(((Xbf.*Kbf)-(Xf.*Kfb)).*S2)))./d; 

CMl=CMa.*.017453; 

CM3=CMa.*.052360; 

CM5=CMa.*.087227; 

CM10=CMa.*.174533; 

CMd=(CNaf.*((Xbf.*kbf)+(Xf*kfb))*S2)./d;-%From pg 20-2 
q=.5.*rho.*((a.*M). x 2); 

CN=pwl./(q.*pcsi) ; 

atrim=((CN.*CMd)./((CNa.*CMd)+(CNd.*CMa))).*57.3; 

LFS=(CNd).*q.*pcsi./pwl; 

cpn=Xn./d; 

cpf=Xf./d; 


figure(1),plot(M,CNa,*k’),grid 

title('CN alpha'),xlabel('Mach'),ylabel(’CN alpha (/rad)') 
figure(2),plot(M,CNd,'k'), grid 

title('CN delta'),xlabel('Mach'),ylabel('CN delta (/rad)') 

figure(3),plot(M,CMa,'k'), grid 

title('CM alpha'),xlabel('Mach'),ylabel('CM alpha (/rad)') 
figure(4),plot(M,CMd,'k'), grid 

title(’CM delta'),xlabel('Mach'),ylabel(’CM delta (/rad)’) 
figure(5),plot(M,atrim,'k') , grid 

title('Alpha trim'),xlabel('Mach'),ylabel('Alpha trim (deg)') 
figure (6).,plot (M,LFS, ' k') ,grid 
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title(* Load Factor Sensitivity 1 ),xlabel(* Mach f ),ylabel( 1 (n-1)/delta T ) 
figure(7),plot(M,CNan, 1 k 1 ),grid 

title{’Cnalpha Nose/Body Alone T ),xlabel(’Mach’),ylabel( 1 CNalpha 
(/rad) 1 ) 

figure(8),plot(M,CNaf,’k 1 ),grid 

title( 1 Cnalpha Fin Alone’),xlabel(’Mach 1 ),ylabel( 1 CNalpha (/rad) 1 ) 
figure(9),plot(M,cpn, 1 k 1 ),grid 

title(’Xcp Nose/Body 1 ),xlabel(’Mach*),ylabel(’Xcp(from eg)/body 
diameter’) 

figure(10),plot(M,cpf, 1 k T ),grid 

title(’Xcp Fin’),xlabel( 1 Mach’),ylabel( 1 Xcp(from eg)/body diameter’) 
figure(11),plot(MjCMl, 1 k’),grid 

title(’CM at One Degree*),xlabel(’Mach’),ylabel(’CM’) 
figure(12),plot(M,CM3,’k’) f grid 

title(’CM at Three Degrees’),xlabel(’Mach’),ylabel(’CM’) 
figure(13) , plot(M,CM5,'k *),grid 

title(’CM at Five Degrees’),xlabel(’Mach’),ylabel(’CM’) 
figure(14),plot(M,CM10,’k’),grid 

title(’CM at Ten Degrees’),xlabel{’Mach’),ylabel( 1 CM’) 


74 



APPENDIX (D) Graph Approximation Programs 


The first program provided approximates the four interference factors. 

% Interference Factors 


clear 

rs=linspace(0,1,20); 


%KBW 
A=.869; 

B=.65; 

C=-.367; 

KBW=(A.*((rs+B). A 2))+C; 

%kBW 

kBW=rs; 

%KWB 

A=.278; 

B=1.3; 

C=.529; 

KWB=(A.*((rs+B). A 2))+C; 


%kWB 
A=. 2; 
B=-.5; 
C=.95; 


kWB=(A.*((rs+B). A 2))+C; 
figure(1) 

plot(rs,KBW,'k',.2,.279, ’ok', .4,.607, ’ok',.6,1,'ok', 
rs,kBW,'k',.2,.213,'dk', .4,.41, 'dk', .6,.61,'dk',. 
rs,KWB,'k',.2,1.16, , *k' / .4,1.35, '*k', .6,1.55, '*k' 


.8,1.446,’ok’, ... 

8,.8,'dk', ... 

,.8,1.77,’*k'),grid 


hold 

plot(rs,kWB,'k',.2,.958, ,A k’,.4,.933, ,A k' 

title(’Interference Factors') 

xlabel('Radius to Semi span Ratio') 

text(.43,1.5,'KWB') 

text(.64,1.3,'KBW') 

text(.21,.91,'kWB') 

text(.41,.35,'kBW') 


6, .94, ’ A k\ .8, .967, ' A k') 


Figure(19) on the next page is the plot of this program 
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Figure (19) Interference Factors 
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The nose/body lift coefficient is calculated by two programs. The curve changes 
when the the fineness ratio equals p and the two branches needed to be done seperately. 

% Nose-Body CNalpha Graphic Approximation 

% Left side of USAF DATCOM figure 4.2.1.1-5a on pg 13-4 of AA3701 notes 

A=-.28; 

B=-l.9; 

C=3.58; 

Beta=linspace(.4,1,100); 

CNalpha=(A.*((Beta+B)- A 2))+C; 
figure(1) 

plot(Beta,CNalpha,'k’,.4,2.95,•ko',.6,3.1,’ko',.8,3.25,'ko',1,3.35, 'ko 
),grid 

title('CNalpha vs Beta/fN') 

xlabel('Beta divided by Nose Length/Body Diameter ratio') 
ylabel('Cnalpha') 


% Nose-Body CNalpha Graphic Approximation 

% Right side of USAF DATCOM figure 4.2.1.1-5a on pg 13-4 of AA3701 
notes 

A=-3.4; 

B=-.9; 

C=3.4; 

Beta=linspace(.1118,1,100) ; 

CNalpha=(A.*((Beta+B). A 2))+C; 
figure(1) 

plot(Beta,CNalpha,'k',.5,2.86,'ko',.6,3.09,'ko',.8,3.36,'ko',.9,3.4, ko 

I 

1,3.35,'ko'),grid 
title{'CNalpha vs Beta/fN') 

xlabel('Beta divided by Nose Length/Body Diameter ratio') 
ylabel('Cnalpha') 


Figure (20) and Figure (21) on the next two pages are the plots of these two programs. 
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Cnalpha 
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The boattail drag and base drag coefficients are calculated from the following 


graphs. 

A=-2000; 

B=65000; 

C=.69; 

D=-570; 

E=l; 

F=.0115; 

G=-.00097; 

M=linspace(1.2, 9,100);' 

m=((D.*M). A 2)-(4.*B.*(<C. *M) + (F. * (M. A 3))+(G.*(M. A 4))+E)); 
cd=(10.*(-(D.*M)-sqrt(m))./(2*B))-.015; 
figure( 1 ) 

plot(M,cd,'k',1.5,.0171,’ ko’ , 2, . 01, ’ ko' , 2.5,.0071,'ko', 3,. 0057, ' ko’, .. 

3.5,.0054,'ko',4,.0051,'ko'),grid 
title(’Boattail Drag Coefficient vs Mach Number’) 
xlabel(’Mach Number’),ylabel(’Boattail Drag Coefficient’) 


% Base Pressure Coefficient Graphic Approximation 
% From Chin pg 31 
A=.007; 

B=1000; 

C=.00007; 

D=-445; 

E=.85; 

F=-.0009; • • 

M=linspace(1,9,100); 

m=((D.*M). A 2)-(4.*B.*((A.*(M. A 2)) + (C.*M)+(F.*(M. A 3)) +E) ) ; 
cd=110.*(-(D.*M)-sqrt(m))./(2*B); 
figure(1) 

plot(M,cd,’k’,1.2,.168,’ko’,1.4,.148, ’ko’,1.8,.12,’ko’,2.2,.1,’ko’,2.6 
.082, ... 

. ’ko’,3.0, .072,’ko’),grid 

title(’Base Pressure Coefficient vs Mach Number’) 
xlabel(’Mach'Number'),ylabel('Base Pressure Coefficient') 


The plots are presented as Figure (22) and Figure (23) on the next two pages. 
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Boattail Drag Coefficient 
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